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The Correlation of Wind-Tunnel and Flight- 
Test Stability and Control Data for an 
SB2C-1 Airplane 


MARION T. HOCKMAN* ann ROBERT E. EISIMINGER, JR.t+ 
Curtiss-Wright Corporation 


SUMMARY 


Adequate correlation of comparable wind-tunnel and flight- 
test data has been neglected until recently. As a result of this 
neglect, much of the accelerated design work that was carried on 
during the war was done purely on the strength of wind-tunnel 
data alone, with little reference to the possibilities of deviation 
from the wind-tunnel results by the actual flying product. In 
some cases, possibly, the differences were small and could be 
neglected, but, in general, extensive programs of revision and 
modification showed up in the aircraft industry on the final flight 
product as a result of failure to correlate the available wind-tunnel 
and flight-test data on earlier models. 

The correlation of comparable wind-tunnel and flight-test 
stability and control data for an SB2C-1 airplane shows the fol- 
lowing results: 


(1) The longitudinal stability determined from flight tests is 
generally greater than that determined from wind-tunnel tests.’ 

(2) The longitudinal control characteristics show fair agree- 
ment as to trends, although there is nearly constant shift in the 
deflection versus velocity curves. 

(3) The stick forces are also in fair agreement, but this is re- 
garded as coincidental, since the effects of the high frictional 
forces found in the airplane control system were neglected in the 
analysis of the wind-tunnel data. In each case the forces were 
found to be of small magnitude throughout most of the flight 
range. 

(4) The elevator effectiveness as determined from the flight- 
test data is slightly greater than that determined from the wind- 
tunnel tests. 

(5) The accelerated flight data show close correlation. 

(6) The correlation obtained for rudder control is not com- 
pletely satisfactory because of the shift in the curves and the 
sometimes large differences in ds,/dV. 

(7) The agreement between the pedal forces is fair, but again 
this is considered coincidental because of the high frictional foices 
found in the airplane control system. 

The source of the discrepancies between the two types of data 
cannot be precisely determined. However, it is believed that the 
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contributing causes can be divided into two groups, each appli- 
cable to one type of data. These are: 


Wind-Tunnel Data 


(1) Scale effects due to the low Reynolds Number of the tests. 

(2) Corrections to the data, such as tare and alignment cor- 
rections or the aerodynamic wall corrections. 

(3) Incorrect duplication of all power effects, such as propeller 
side forces, because of the use of a fixed-pitch propeller on the 
model. 

(4) Omission of some external drag items on the model, such 
as landing gear, lap joints, etc. 

(5) Use of a solid wooden model having lower percentage 
structural deflections under load than the airplane. 


Flight-Test Data 


(1) Incorrect estimation of the c.g. position, which must, of 
necessity, be based on an average c.g. position because of fuel 
consumption during any flight. 

(2) Presence of any longitudinal or vertical accelerations 
when the data are recorded. 

(3) Variation in flying technique of different pilots, or day-to- 
day variation in an individual pilot. 

(4) Deflection of airplane structures under air loads during 
flight, particularly wing or tail surfaces. 

(5) Variations.in instrument accuracy and calibration due to 
atmospheric conditions, accelerations, or vibration. 


It is not possible to state which of these items have the great- 
est effect, since some will undoubtedly be of secondary impor- 
tance. In view of the many possible sources of error, it is con- 
sidered that the results of the correlation presented in this analy- 
sis are generally satisfactory and will further establish the use of 
wind-tunnel tests to predict the flight characteristics of airplanes. 


NOMENCLATURE 


airplane weight, Ibs. 

area of wing, sq.ft. 

area of aileron, sq.ft. 

area of elevator, sq.ft. 

‘area of rudder, sq.ft. 

wing mean aerodynamic chord (M.A.C.), in. 
center of gravity movement, in. 
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Ce = aileron mean geometric chord aft aileron hinge line, TABLE I 
»  $B2C-1—Physical Characteristics 
Ss = elevator mean geometric chord aft elevator hinge 
49.719 ft. 
or = rudder mean geometric chord aft rudder hinge line, y= (projected in root chord plane) 422 sq.ft. 
2 in. Area bla ted by fuselage 57.9 sq.ft. 
Root chord (at centerline of airplane) 144 in. 
B = angle of sideslip, positive when right wing is for- Construction tip chord 62 in, 
ward, deg. ent tip chord 63. in 
n aerodynamic chord 109.3 iu. 
= angle of yaw itive w ing is forwar of airplane to M.A.C. (lateral) 125.56 in, 
v & yaw, positive hen left wing 1 f d, Root chord plane to M.A.C. (vertical) 13. 12 ia 
deg. of wing to chord point 
26.67 
v = og speed, ft. per sec. Distance from leading edge of wing aft to ™ 
Vm.p.b. = air speed, miles per hour 
p = mass air density, slugs per cu.ft. Incidence—root +1. 3 
o = density ratio pay 
nm span é 

q = dynamic pressure, lbs. per sq.ft. [¢ = (?/2p)v?] Average aileron chord, aft of hinge line 17.63 in 
C if ffici - Aileron area (total), aft of hinge line 26.6 = sq.ft. 
L = lift coefficient (lift /gS) incl. balance 23.48 in. 

Cn = pitching — coefficient about c.g. (pitching Aileron balance area 3'8 ae 
momen 4 ea affecte vy aileron 120 ft. 
q Airfoil section—root_- NACA 23017 . 
Ch = hinge moment coefficient. Subscript a, e, or r de- Airfoil section—tip NACA 23009 
notes aileron, elevator, or rudder. Based on sat 
mean geometric chord and control surface area _— affected by upper flap 205° sq.ft 
F, = stick force, Ibe. by lower flap sq.ft. 
Sweepback (leading edge) 
Fy pedal Ibs. Slats Partial Span Handley-Page 
é = control su ace setting, deg. Subscript a, e, or r (B) Horisontal Tail Surfaces 
denotes aileron, elevator, or rudder Total Area 107.4 sq.ft. 
Span 228.5 in. 
As = total angular range of control surface, deg. Sub-— 9157 
script a, e, or r denotes aileron, elevator, or rud- Area affected by elevators 95.25 sq.ft 
d Maximum chord 99.84 in. 
er Stabilizer area (forward of pen H.L.) 79.6 sq.ft. 
Ad = angular stick or pedal displacement corresponding in 
to Aé, deg. (aft of H.L.) 18.5 in. 
evator balance sq.ft. 
= sti-k length, in. Tab area (total) 3 2.84 
lp = pedal arm length, in. Modified N-69 
dC,,/d;, = change in pitching moment with change in stabil- (C) Vertical Tail Surfaces 
izer incidence Total area 45.7  sq.f 
5. q.ft. 
= change in elevator hinge moment with change in 60:0 safe 
stabilizer incidence (to 26.5 sq.ft 
ds,/dy = change in rudder deflection required for trim per Rudder height 10712 
verage rudder chord (aft o 25.8 in. 
degree change in yaw angle 14 4 sqft. 
n = load factor ).85 sq.ft. 
= acceleration of gravity, ft. per sec. per sec. or 
£ = tunnel wall correction factor, equal to 0.105 for (p) Engine Data 
Wright Brothers Wind Tunnel Mfg. design R-2600-8 - 
= down-wash factor for tunnel wall correction to 
Supercharger 2 Speed 
pitching moment Ratings 1 700, hp. at 2,800 r.p.m. at 
A = tunnel cross-section area, 57.3 sq.ft. 1, 00" hp. at 2,400 r.p.m. at 
at 2,400 r.p.m., at 
INTRODUCTION 350 bp. at 2,400 r.p.m, at 


La wade BY THE IMPETUs of urgent wartime require- 
ments, many combat airplanes were designed and 
put in production before adequate flight-test evaluation 
of the handling qualities could be completed. Conse- 
quently, the designer was forced to accept the results 
of wind-tunnel tests as the basis for any changes re- 
quired to ensure satisfactory qualities. Since the pilot 
becomes the final judge of the airplane, it is and should 
be the goal of the designer to supply him, in the shortest 
possible time, with an article that he considers satis- 
factory, particularly in regard to the handling qualities. 
Unfortunately, there is considerable disagreement 
among pilots as to what constitutes a satisfactory air- 
plane, and the designer must rely on a general stability 
and control specification such as that developed by the 
Navy in reference 8. 

Obviously, if the designer is forced to rely on the 
results of wind-tunnel tests to justify any design 


(E) Propeller Data 

Type Curtiss 
Blade Design 89324-12 
Diameter 12 ft. 
Number of blades 


(F) Control Movement 


Total movement of stick for elevators 40° or E = 55°) 


Length of stick for elevators 231/2 i 
Total movement of rudder pedals 31° 33° ‘(for R = 50°) 
Radius of pedal movement about center of 
rotation 113/, in 
Total movement of stick for ailerons 0° (for A = 25°) 
«in 


Length of stick for ailerons 


changes, it is necessary to know how reliable are the 
data on which any analysis may be based. The first 
step seems to be a correlation of comparable wind- 
tunnel and flight-test data. 

In order to better determine the reliability of powered 
wind-tunnel tests in predicting stability and control 
characteristics of airplanes, some comparable wind- 
tunnel and flight-test data for a Curtiss SB2C-1 airplane 
were analyzed so that the degree of correlation between 
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STABILITY AND CONTROL DATA FOR THE SB2C-1 


the two types of data could be readily ascertained. An 
extensive instrumented flight-test stability and control 
investigation was conducted by the N.A.C.A. on an 
SB2C-1 airplane, Serial No. 00014, at Langley Field, 
Va., and the results are presented in reference 1. 

After the flight tests had been made, wind-tunnel 
tests were conducted on a one-eighth scale model of 
the SB2C-1 airplane at the Wright Brothers Wind 
Tunnel at the Massachusetts Institute of Technology, 
as part of a general program studying the effects of 
various modifications on the size, weight, performance, 
and flying qualities of the SB2C-1 airplane. For this 
reason the wind-tunnel tests did not completely dupli- 
cate the flight-test program. Consequently, it will 
not be possible to correlate some of the lateral and direc- 
tional control data’ The results of the wind-tunnel 
tests are given in reference 2. This analysis should be 
useful because it indicates the accuracy and reliability 
of the tunnel data in the estimation of the flying quali- 
ties for this particular airplane, although it must be 
emphasized that any discrepancies that may be indi- 
cated are not intended as a criticism of the Wright 
Brothers Wind Tunnel. It is probable that different 
results would be obtained from other data. It is in- 
tended to point out the desirability of similar correla- 
tions for other airplanes whenever comparable data are 
available. 


DESCRIPTION 


The SB2C-1 airplane is a two-place, single-engined, 
low-wing cantilever monoplane with a retractable land- 
ing gear, partial span split flaps, and Handley-Page tip 
slats. The flaps are double split, so that they may be 
used either as landing flaps or dive brakes. The slats 
are located along the leading edge of the wing in front 
of the ailerons and are extended automatically with the 
landing gear. The airplane is equipped with the Wright 
R-2600-8 engine and a three-bladed Curtiss constant- 
speed propeller of 12.17 ft. diametet. The control sur- 
faces are fabric covered. The stabilizer incidence is 
+3° with respect to the thrust line, and the fin setting 
is —1.5° with respect to the airplane centerline. Addi- 
tional characteristics can be found in Table 1. The 
general airplane arrangement is shown in Fig. 1. 

The one-eighth scale model was of wood construction 
and was equipped with an electric synchronous motor 
and propeller for powered testing. The surface of the 
model was finished with pigmented lacquer and rubbed 
down to a smooth surface. An orifice plate was 
mounted in the cowling, designed for a pressure drop 
that was intended to duplicate the flow characteristics 
through the engine compartment only for the high- 
speed condition. Many of the detailed items on the air- 
plane, such as antenna masts, lapped skin junctions, 
landing gear for the landing configuration, etc., were 
omitted on the model. The model propeller was a fixed 
pitch unit. The model stabilizer incidence was +3° to 


Fic. 1. General airplane arrangement. 


the thrust line, and the fin incidence was —1.5° with 
the model centerline. 


CORRECTIONS TO DaTA 


The standard Glauert wall corrections to drag and 
angle of attack were applied to the wind-tunnel data. 

Pitching moment corrections to the wind-tunnel 
data due to wall restraint on the flow at the tail were 
computed from 

S dC, 

AC,, = —57.3§r r (1) 
where the factor 7 is a function of the tail length and 
was derived by Lotz.’ 

Tare and alignment corrections were incorporated in 
the wind-tunnel data as given in reference 2. No cor- 
rections were applied to the wind-tunnel control surface 
hinge moments., 

The flight-test instruments were calibrated on the 
ground, and the necessary corrections were applied to 
the flight-test data. Where necessary, position error 
corrections were also applied. 


PROCEDURES 


General 


In analyzing the data, a gross weight of 12,389 Ibs. 
was used. Reference 1 did not indicate the gross weight 
for each c.g. position, and, therefore, the above gross 
weight is an average of the gross weight range of 12,000 
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Fic. 2. Pitching moment transfer for horizontal c.g. movement. 
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Fic. 3. Relation between control surface deflection and stick and 
rudder position. 


to 12,672 lbs., with a corresponding c.g. range of 23.3 to 
32.3 per cent M.A.C., used during the flight-test pro- 
gram. When referring to the data, indicated velocities 
will be used throughout the report. Some error, not 
serious however, is introduced in this analysis, since the 
wind-tunnel power data were based on a gross weight of 
13,509 Ibs. and sea-level altitude. 

The flight-test and wind-tunnel control forces were 
trimmed to the velocities indicated in the following 
table: 
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Configuration Velocity 
Cruising, power-off 300 m.p-h. 
Cruising, rated power 300 m.p.h. 
Landing, power-off ‘96 m.p.h. 
Landing, rated power 100 m.p.h. 
Diving, power-off 300 m.p.h. 


The methods used in trimming the data will be ex- 
plained later. 


Analysis of Wind-Tunnel Data 


The indicated velocity corresponding to a given lift 
coefficient was computed from 


V, = V391W/SC, (2) 


The control surface deflections were measured with 
respect to the stabilizer and fin chord line, respectively, 
during the test program. 

The stick fixed neutral point variation throughout 
the flight range was computed from the stabilizer ef- 
fectiveness data. The graphical method given in ref- 
erence 4 was used for this computation. 

The wind-tunnel pitching moment curves were plot- 
ted directly from the tabulated data of reference 2, 
from which it was necessary to transfer the data to the 
flight-test c.g. positions. This was done approximately 
by the following method. Referring to Fig. 2, the equa- 
tion for the moment transfer from (c.g.); to (c.g.)2 is 


M, = M, + Lacos a + Dasin a (3) 


or in coefficient form 
a a 
Cm + G, cos a + Cy sin a (4) 


By approximation, cos a = 1 and sin a = 0; therefore 
Cm = Cm + (a/c)Cy (5) 


Eq. (5) was used for transferring the wind-tunnel 
pitching moment data to the flight-test c.g. positions. 

The longitudinal control characteristics were ob- 
tained from the elevator effectiveness data in reference 
2. The stick forces were computed from the wind- 
tunnel data using the following control force equa- 
tion: 


le Aé, 1 (6) 


Since there were no wind-tunnel tab effectiveness data, 
the flight-test tab effectiveness was used in determining 
the required elevator trim tab angle at the specified 
trim velocity. 

The stick force and elevator deflection required for 
the accelerated flight condition were computed by the 
method of reference 5. 

The directional control characteristics were deter- 
mined from the wind-tunnel rudder effectiveness data 
given in reference 2. The pedal forces were calculated 
using the following equation: 
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Aé, 1 
— (7) 


F, =— 


The flight-test tab effectiveness was used in determining 
the required rudder trim tab angle at the specified trim 
velocity. 

In the control force equations, the values of the me- 
chanical advantage factor, Ad/A@, were assumed to re- 
main constant, irrespective of control surface deflection. 
That this assumption is reasonable can be seen in Fig. 
3, which shows that the stick and rudder positions 
closely follow a linear relationship with the control sur- 
face deflections. The data for Fig. 3 were obtained by 
actual measurement on the ground. 

The wind-tunnel directional control data were modi- 
fied to include the flight-test angle of sideslip shown in 
Figs. 19 through 22. This was done by determining a 
value for the parameter, d6,/dy, from the wind-tunnel 
tests. The angle of yaw is equal in magnitude to the 
angle of sideslip but opposite in sign. Therefore, 


from which the incremental change in rudder angle re- 
quired to trim when the airplane is flying at small 
sideslip angles is 


Aé, = —(di,/dp) AB (9) 


Applying this increment to the wind-tunnel data places 
the wind-tunnel directional control data on a compa- 
rable basis with the flight-test data. 

The directional control data were not transferred 
to c.g. positions corresponding to those used during the 
flight tests, since at or near zero sideslip any such cor- 
rections would be small. Neither were the data cor- 
rected for the yawing moment induced by the deflected 
aileron required to hold the wings level. , An examina- 
tion of the data of references 1 and 2 indicated this cor- 
rection to be of negligible magnitude, since the aileron 
deflection required always remained small. 


Analysis of Flight-Test Data 


For the flight tests, reference 1, the control position 
indicators were rigged so that the elevator angles were 
measured with respect to the thrust line and the rudder 
angles were measured with respect to the airplane cen- 
terline. The deflections for the elevator and rudder 
have been corrected to measure from the stabilizer and 
fin chord lines, respectively, for this analysis, since this 
is the normal production rigging. 

The lift coefficient corresponding to a given velocity 
was computed from 


C, = 391W/SV? (10) 


The stick fixed neutral point variations were com- 
puted from the flight-test data by the use of the eleva- 
tor control data. The method used for determining the 
stick fixed neutral point variation is given in reference 
6. To apply the method, a plot of elevator angle, Se, 
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Fic. 4. Variation of elevator deflection with lift coefficient, 
cruising configuration, power-off. 
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Fic. 5. Variation of elevator deflection with lift coefficient, 
cruising configuration, rated power. 


versus the airplane lift coefficient, C,, for at least three 


center of gravity positions is required; Figs. 4 and 5 
are typical and were copied from reference 1. 

The airplane pitching moments are obtained from 
the flight-test longitudina! control data by the following 
method. The method consists of determining the C;’s 
for trim from the several c.g.’s tested (at least three 
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c.g.’s are necessary) at a constant elevator angle. 
These points are spotted on the C, scale for C,, = Oofa 
Cm versus C, graph. The change in C,, due to a change 
in c.g. from c.g.; to c.g.2 is 


— 
(11) 


A curve of C, vs. Cm, can now be constructed as shown 
in Figs. 11 and 12. 

When pitching moments are desired for elevator 
angles for which there are no trim lift coefficients at any 
of the several c.g.’s tested, it is necessary to obtain 
values for the parameters dC,,/d6, for the flight range. 

For a constant C;, on the 6, versus C, plot, Fig. 5, a 
change in c.g. causes a change in the elevator deflec- 
tion required for trim. Therefore, a AC,, can be com- 


puted from 


AC, = we Cm = — 


AC, = —(Ac.g./100) Cz, (12) 
The corresponding Aé, can also be computed 
Aé, = (13) 


Dividing Eq. (12) by Eq. (13) gives the dC,,/dé, for 
the given C,, 


Ac.g. 


This is done for several values of lift coefficient in order 
to determine the elevator effectiveness throughout the 
flight range. Having previously computed a pitching 
moment curve for one elevator setting, usually zero 
degrees, additional curves can be calculated for the de- 
sired elevator deflection. 

This parameter will be satisfactory for low elevator 
deflections, say +10°. For elevator deflections beyond 
this range, the value of the parameter may decrease 
because the elevator effectiveness tends to decrease. 

The flight-test and wind-tunnel longitudinal control 
data, Figs. 14 through 18, were trimmed at the veloci- 
ties indicated previously. This was accomplished by ob- 
taining the elevator tab effectiveness at the required 
trim velocity from Fig. 13. The amount of tab required 
was determined as follows: 

At the required trim velocity, obtain dF,/dé, from 
Fig. 13. Then, 


(Fs) ser0 (15) 


The stick force at the trim velocity will now be zero. 
A new force curve can then be obtained for the fixed 
-trim tab setting throughout the velocity range by the 
use of Fig. 13. 

The maneuvering flight data were plotted directly 
from reference 1. The major portions of the data were 
taken at an air speed of 200 m.p.h. However, there 
seemed to be a negligible change in the data between 200 
and 280 m.p.h. It is estimated that there would be no ap- 


preciable difference between 200 and 240 m.p.h. An Fic.9. Neutral point analysis, landing configuration, rated power. 
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Fic. 6. Neutral point analysis, cruising configuration, power-off. 
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Fic. 7. Neutral point analysis, cruising configuration, rated 
power. 
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Fic. 8. Neutral point analysis, landing configuration, power-off. 
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air speed of 240 m.p.h. is used because the wind-tunnel 
data were computed for this air speed. 

The flight-test and wind-tunnel directional control 
data were trimmed by the same method as that used in 
trimming the longitudinal control data. 


DISCUSSION 


The similarity of the configuration as tested in the 
wind tunnel and in flight must be compared before any 
attempt can be made to coordinate the results between 
the two types of tests. Table 1 will clarify the reference 
to the various configurations later in the discussion. 


Longitudinal Stability and Control 


The neutral point comparison indicates that the 
wind-tunnel data are generally conservative. For the 
cruising configuration with power-off, Fig. 6, the flight- 
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Fic. 10. Neutral point analysis, diving configuration, power-off. 
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Fic. 12. Variation of pitching moment coefficient with lift 
coefficient, cruising configuration, rated power. C.g. = 24 
per cent M.A.C. Tail incidence = 3°. 


test data show a 3'/2 per cent stability margin over the 
wind-tunnel data in the Ver, and Vain,» range. In 
the Vyiae range, the two curves converge, and the wind- 
tunnel data indicate better stability than the flight- 

test data. 

The application of rated power for the cruising con- 
figuration, Fig. 7, is destabilizing in each case, but the 
flight-test data maintains a 1 per cent margin of sta- 
bility over the wind-tunnel data throughout most of 
the speed range. For high speed the flight-test and 
wind-tunnel neutral point locations are the same. 

For the landing configuration with power-off, Fig. 8, 
the flight-test data show a 3 per cent average stability 
margin over the wind-tunnel data throughout the flight 
range. Rated power, Fig. 9, tends to bring the flight- 
test and wind-tunnel data in close agreement, except 
for high approach velocities where the wind-tunnel 
data indicate better stability. 

For the diving configuration, Fig. 10, the airplane 
shows a 3 per cent margin in stability over the model. 
However, the stability trend diverges sharply at lift 
coefficients greater than one, the reason for which is 
not readily apparent. 

In connection with the neutral point data it was con- 
sidered advisable to compare the pitching moments as 
obtained from the flight-test and wind-tunnel data. 
The pitching moment curves, Figs. 11 and 12, show the 
same degree of comparison as the neutral point analy- 
sis, which is to be expected since both were based on the 
same data. However, it is to be noted that the elevator 
effectiveness determined from the flight-test data is 
slightly greater than that shown by the wind-tunnel 
data. 


a 
| | | 
| 
i 
% 

| 

off. 

| 
| 
| | 
a 
| 
i 
x | 
" : 
2.2 


: 12 JOURNAL OF THE AERONAUTICAL SCIENCES—JANUARY, 1948 
zw 
CRUISING CONFIGURATION 
<0 
POWER OFF INDICATED VELOCITY, MPH 
60 100. 40 300 340 20 
INDICATED VELOCITY, M. H. 
+ 
i RUDDER TAB 
| | | + | | | 
wa Re | (FLIGHT TEST RESULTS) 20 WIND TUNNEL#2.31°TAIL HEAVY | 
100-40 300 340 | | | 
INDICATED VELOCITY,M.P.H. 
Fic. 13. Control force variation with trim tab angle. co 
Fic. 15. Elevator deflection and stick force for trim, cruising ” 
configuration, rated power. C.g. = 24.9 percent M.A.C. Trim 
velocity = 300 m.p.h. 
go 
20 da 
| tic 
| ity 
fec 
8 cu 
t- CE WIND TUNNEL RESULTS | 
WwW + - 4 
8 \ | | T = 
WIND TUNNEL RESULTS 80 100 120 140 160 180 
100 140 180 340 ELEVATOR TAB ANGLES FOR TRIM 
IMOIGATED VELOCITY, MPL. 20 | FLIGHT TEST= 20.3°TAIL HEAVY - 
rat | WIND TUNNEL= 3.2°TAIL HEAVY 
ELEVATOR TAB ANGLES FOR TRIM 
ol FLIGHT TEST = 1.4° TAIL HEAVY — re) 
WIND TUNNEL = TAIL HEAVY 
| WIND TUNNEL RESULTS So 
” a | | 
10, | ~ 
INDIGATED VELOGITY,. M.P.H. 
Fic. 14. Elevator deflection and stick force for trim, cruising Fic. 16. Elevator deflection and stick force for trim, landing 
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Fic. 17. Elevator deflection and stick force for trim, landing 
configuration, rated power. C.g. = 24.0 percent M.A.C. Trim 
velocity = 100 m.p.h. 


The longitudinal control data, Figs. 14 and 18, show 
good agreement for all configurations. However, for the 
cruising configuration with power-off, the flight-test 
data show a slightly greater change in elevator deflec- 
tion with velocity, indicating that the increased stabil- 
ity has a greater effect than the increased elevator ef- 
fectiveness. With rated power the elevator deflection 
curves, though displaced, are nearly parallel, thereby 
indicating that the increased stability and elevator ef- 
fectiveness of the flight-test data tend to cancel each 
other. 

The stick forces likewiseshowfair agreement, although 
this is regarded as more or less coincidental because of 
the 5-lb. friction force present in the airplane control 
system, for which no allowance was made in the analysis 
of the wind-tunnel data. In many cases the friction 
force was greater than the measured stick force. This 
may partially explain the differences in trim tab angles 
required. 

There were no tab effectiveness data obtained during 
the wind-tunnel tests, and it was necessary to use the 
flight-test tab effectiveness data from Fig. 13 in deter- 
mining the tab angle for trim for the wind-tunnel data. 
The use of the flight-test tab effectiveness in the analy- 
sis of the wind-tunnel data is questionable, since it is 
possible that the wind-tunnel tab effectiveness would 
be different from that measured in flight test. However, 
it is considered more desirable than the use of calcu- 
lated tab effectiveness. 

A comparison of the maneuvering flight data is shown 
in Figs. 19 and 20. The variation of control position 
with lift coefficient and stick force with load factor is 
shown in Fig. 19 for a c.g. position of 28 per cent 
M.A.C. The maneuvering neutral point, stick-fixed 
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Fic. 18. Elevator deflection and stick force for trim, diving 
configuration, powar-<e. C. g. = 24.4 per cent M.A.C. Trim 
velocity = 300 m.p.h 


12, 
uj 8 
4|_WIND TUNNEL RESULTS 
r 
— LIGHT TEST RESULTS 
% 4 Ke) 
LIFT COEFFICIENT, CL 
3 
A 
Pio TUN 
x | RESULTS 
b | | 
% 1 4 5 


LOAD FACTOR, n 


Fig. 19. Longitudinal characteristics in steady turns, cruis- 
ing configuration, power-off. C.G. = 28 per cent M.A.C. 
Trim velocity = 240 m.p.h. 


and stick-free, is shown in Fig. 20. These data indicate 
that the correlation is well within satisfactory limits 
based on the magnitude of control system friction that 
was inherent in the flight airplane. 

It may be of interest to compare the hinge moment 
parameters dC,,/da, and dC,,/dé, as obtained from the 
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Fic. 20. Maneuvering neutral point, cruising configuration, 
power-off. Trim velocity = 240 m.p.h. 
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Fic. 21. Rudder deflection and pedal force for ttim, cruising 
configuration, power-off. C.g. = 24.55 per cent M.A.C. Trim 
velocity = 300 m.p.h. 


flight and wind-tunnel tests. The values of these pa- 
rameters were determined to be 


Wind Tunnel Flight Test 
0 —0.0012 
—0.0024 —0.0033 


The comparison is not, in general, considered com- 
pletely satisfactory. 


Directional Control 


The comparison between directional control data for 
the cruising configuration, Figs. 21 and 22, is satis- 
factory. The cruising configuration, power-off, indi- 
cates a difference in rudder deflection required for 
trim of 1.3° throughout the speed range. However, 
the control gradient, dé,/dV, is approximately the same. 
Rated power causes a steeper control gradient for the 
wind-tunnel data than for the flight-test data. 


There is a discrepancy of considerable magnitude in 
the rudder required for trim for the landing configura- 
tion, Figs. 23 and 24. With power off the difference 
varies from 2'/2° at 100 m.p.h. to 5'/2° at 160 m.p.h. 
The difference using rated power is approximately 
5'/2° throughout the flight range. If the flight-test 
data had not been available prior to the wind-tunnel 
tests, it is possible that considerable effort would have 
been expended in an attempt to improve the directional 
control in a wave off, which is characterized by high 
powers at very low speed and is generally the most 
critical from a control standpoint. For this reason, the 
wind-tunnel data for this condition is not considered 
satisfactory. 


The directional control forces do not agree in magni- 
tude at various air speeds. However, the trends are 
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Fic. 22. Rudder deflections and pedal force for trim, cruising 
configuration, rated power. C.g. = 24.9 per cent M.A.C. Trim 
velocity = 300 m.p.h. 
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similar for the wind-tunnel and flight-test data through- 
out the flight range for all configurations. Again this is 
regarded as coincidental because of the large friction 
forces found on the airplane. The friction force varia- 
tion with rudder deflection is shown in Fig. 25. 

The rudder tab angles required for trim are subject 
to the same discrepancies as the elevator tab angles, 
since the flight-test tab effectiveness was used to deter- 
mine the tab angle for trim for the wind-tunnel data. 

There were no comparable directional control data 
available from the wind-tunnel tests for the diving 
configuration. However, it is believed that a compari- 
son of data for the diving configuration would agree as 
well as that shown for the cruising configuration, since 
the airplane flight characteristics are similar for these 
two configurations. 

The values of the rudder hinge moment parameters as 
determined from the wind-tunnel and flight-test data 
are: 

Wind Tunnel 


doe 0 
dCn,/d6, —0.0028 


Flight Test 
0 
—0.0028 


The comparison is in much better agreement than 
that for the elevator. 

No data are shown for directional control at various 
angles of sideslip, since there were not sufficient wind- 
tunnel data available. What data were available did 
not represent the correct lift coefficients. Similarly, 
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Fic. 23. Rudder deflection and pedal force for trim, landing 
configuration, power-off. C.g. = 23.7 per cent M.A.C. Trim 
velocity = 96 m.p.h, 
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Fic. 24. Rudder deflection and pedal force for trim, landing 
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Fic. 25. Variation of rudder friction force with pedal position. 


no comparison is shown for lateral control because of 
the incomplete wind-tunnel data. 


DIFFERENCES IN THE DATA 


The wind-tunnel and flight-test data are subject to 
various differences, some of which may be apparent in 
the data and some which are considered likely to exist 
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TABLE 2 
Comparison of Flight Test and Wind Tunnel Configurations 


(I) Cruising Configuration—Power-Off* 
Wind Tunnel 


Landing flaps retracted 
Slats retracted 
Landing gear retracted 
Dive brakes retracted Dive brakes retracted 
Cowl flaps closed Cowl flaps closed 


(II) Cruising Configuration—Rated Power* 
Wind Tunnel 


Same as the power-off configu- 
ration 


_ Flight Test 
Landing flaps retracted 
Slats retracted 
Landing gear retracted 


Flight Test 
Same as the power-off con- 
figuration except the cowl 
flaps were opent 


(III) Landing Configuration—Power-Off* 
Wind Tunnel 


Landing flaps extended (60°) 
Slats extended 

Landing gear retracted 
Dive brakes retracted 

Cowl flaps closed Cowl flaps closed 
Front cabin closed Front cabin closed . 


(IV) Landing Configuration—Rated Power* 
Wind Tunnel Flight Test 
Same as the power-off con- Same as the power-off con- 
figuration figuration except the cowl 
flaps were opent 


(V) Diving Configuration—Power-Of* 
Wind Tunnel 
Landing flaps retracted 
Slats retracted 
Landing gear retracted 
Dive brakes extended (upper 
and lower 25°) 
Cowl flaps closed 
Front cabin closed t 


* Power-off for the model is a windmilling condition in which 
no power is drawn from the model motor, while for the airplane 
power-off is an idling condition in which the propeller is drawing 
some power from the engine. Rated power for both wind tunnel 
and flight test.correspond to 1,500 hp. at sea level. However, 
rated power at sea level is greater than rated power at altitude, 
for altitudes above the engine critical altitude. 

t+ There was no landing gear or movable cabin on the wind- 
tunnel model. 

¢t There were no cowl flaps on the model. Cowl flaps were 
opened on the airplane in order to obtain proper engine cooling. 


Flight Test 
Landing flaps extended (60°) 
Slats extended 
Landing gear extended 
Dive brakes retracted 


Flight Test 
Landing flaps retracted 
Slats retracted 
Landing gear retracted 
Dive brakes extended (upper 
and lower 25°) 
Cowl flaps closed 
Front cabin open 


but whose magnitudes are unknown. These differences 
can best be explained by listing and briefly discussing 
those applicable to each type of data. 


Wind-Tunnel Data 


(1) Scale Effects——Because the airplane velocity is 
greater than the velocity at which tunnel tests are made, 
the Reynolds Number of the airplane is greater, pro- 
vided the altitude at which the airplane is flown is not 
excessive. The initial flow separation on the model 
wing, occurring at a lower angle of attack than on the 
airplane, should manifest itself as a progressive increase 
in longitudinal stability up to the model stall because 
of the breakdown in down-wash angle at the tail. The 
divergence of the neutral point data for the diving con- 
figuration is unexplainable at this time because there 
were not sufficient data available to determine the 
cause. Principally, scale effects will cause the greatest 
differences in maximum lift coefficient and longitudinal 
stability near the stall, provided that the angle of at- 


JOURNAL OF THE AERONAUTICAL SCIENCES—JANUARY, 1948 


tack of the tail surfaces is not large enough to result in 
a tail stall at lower lift coefficients. 

(2) Corrections to the Data.—The two types of cor- 
rections applied to the data consist of a mechanical cor- 
rection due to flow misalignment and balance tares and 
an aerodynamic correction due to the effect of the tun- 
nel walls on the model characteristics. The mechanical 
correction is generally determined for the propeller-off 
condition and is probably not correct for a power-on 
condition. 

The classical wind-tunnel wall corrections are based 
on an ideal wing having an elliptical lift distribution 
and are not, therefore, strictly applicable to any other 
wing or to a complete model, although it can be shown 
that minor variations from an elliptical lift distribution 
will have a negligible effect on the magnitude of the 
correction. The application of power may cause 
abrupt changes in the lift distribution because of slip- 
stream twist and velocity and should, therefore, require 
an additional correction of the type derived by Swan- 
son.’ 

No wall corrections are applied to rolling or yawing 
moments or control surface hinge moments, since such 
corrections are small and generally negligible. | 

(3) Duplication of Power Effects—The propeller on 
the model is a fixed-pitch unit and does not, therefore, 
duplicate completely the characteristics of a full-scale 
constant-speed propeller, the most important of which 
may be the propeller side forces. Some calculations 
made in conjunction with another model indicated that 
the side-force coefficients of the model propeller were 
nearly three times those of the full-scale propeller. 
However, a calculated correction for this difference can- 
not be justified. 

Similarly, because the model and airplane torque coef- 
ficients are not equal, the slipstream twist angles are 
not equal, which will manifest itself in the yawing 
moment data. Interpretation of the data may then 
give misleading results. 

(4) Omission of External Drag Items on the Model.— 
It is customary to omit from the model numerous exter- 
nal drag items such as rivets, lap joints, radio masts, 
etc., since such items have been found to be extremely 
sensitive to changes in Reynolds Number. This omis- 
sion is satisfactory provided the items are not located on 
the wing where induced flow separation could be started 
which could affect stability characteristics; also, if 
the drag items are large, the stability may be changed 
because of the drag moment. The latter is particularly 
true in this analysis, since the landing gear and cowl 
flaps were omitted from the model. 


Flight-Test Data 


(1) Estimation of Center of Gravity.—Since all of the 
longitudinal data are critically affected by the c.g. posi- 
tion, it is necessary to determine as accurately as pos- 
sible the flight-test c.g. position. Because of continu- 
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ous consumption of the fuel load, this is sometimes dif- 
ficult to do, and it is often necessary to use an average 
c.g. position. For planes whose c.g. is not affected 
by fuel consumption, this problem is of little impor- 
tance. 

The effects of c.g. movement due to fuel consumption 
on directional gontrol data at or near zero yaw are small 
and need not be considered in this analysis. 

(2) Presence of Longitudinal or Vertical Acceleration 
in Data.—To evaluate static stability correctly it is 
desirable that no acceleration other than that due to 
gravity be present in the data. Provided that a precise 
flying technique is carefully followed, this will normally 
be true, but gusts or other unusual weather conditions 
may prevent satisfactory data from being obtained. 
In addition, when steady state longitudinal data are ob- 
tained, it is desirable that the rate of change of velocity 
with time be not greater than approximately '/: m.p.h. 
per sec. 

(3) Variation in Flying Technique.—One of the most 
difficult items to evaluate is the variation in flying 
technique of different pilots or the day-to-day variation 
in an individual pilot, since the accuracy and usefulness 
of the data will depend almost entirely on the pilot. 
For this reason it is advisable, if possible, for one pilot 
to make all the tests or, failing this, to have pilots of at 
least equal abilities and temperament. It is not possible 
to state what influence this may have on the analysis 
in this report. 

(4) Deflection of Airplane Struciures in Flight.— 
The airplane structure, particularly the wings and tail 
surfaces, is elastic enough to cause noticeable changes in 
all stability parameters because of deflections. Because 
the wind-tunnel model is much more rigid, the deflec- 
tions are smaller and the effect on the stability param- 
eters is less. These effects undoubtedly may cause 
some of the discrepancy shown between the wind-tunnel 
and the flight-test data. 

The control surfaces provide a graphic example of 
the effect of air loads on surface deformation, since it 
has been demonstrated that the effects of fabric bulge 
under load may cause as much as 100 per cent variation 
in the hinge moment characteristics of closely balanced 
surfaces. The use of a solid control surface for the 
prediction of stick and pedal forces is therefore of doubt- 
ful value. 

(5) Changes in Instrument  Calibration.—-The 
changes in instrument calibrations or malfunctioning 
of an instrument during any flight is often unknown 
until the data are ready for reduction, at which time it 
is necessary to either neglect the data for that instru- 
ment or to repeat the flight. The wind-tunnel tests are 


more closely supervised, and a running check of the 
data as taken generally reveals any discrepancy. 

There is also some doubt about the accuracy of the 
method for obtaining some of the flight-test data—for 
example, the measurement of the sideslip angles. The 
sideslip angles are determined by means of a vane at- 
tached to a mast projecting ahead of the leading edge of 
the wing at the tip, and it is assumed that the equilib- 
rium angle of the vane is the sideslip angle. Such a 
method probably does not possess sufficient accuracy to 
be used for a close comparison with wind-tunnel data 
because of the influence of the field of flow ahead of the 
wing on the airflow direction around the vane. 


CONCLUSIONS 


The correlation of comparable wind-tunnel and flight- 
test data for an SB2C-1 airplane verifies the reliability 
and usefulness of powered wind-tunnel tests for pre- 
dicting stability and control characteristics of airplanes. 
Several basic differences in the data have been briefly 
discussed, and the conclusion is reached that the degree 
of correlation for this airplane and model is satisfactory. 
However, correlation of the data for one model is not 
considered satisfactory. The correlated data of several 
airplanes might well be examined statistically in order 
to determine to what degree wind-tunnel data may be 
trusted in the design of future airplanes. d 
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Buckling of Curved Sheet in Compression 
and Its Relation to the Secant Modulus 


E. H. SCHUETTE* 


The Dow Chemical Company |. , 


SUMMARY 


Compression tests were made on 87 curved magnesium-alloy 
panels with aspect ratios near unity to determine the buckling 
stress. Panel radius-thickness ratios ranged from 86 to 515. It 
was found that the buckling stress for such panels could be pre- 
dicted with equivalent accuracy within and beyond the elastic 
range by the formula f.,, = 0.42E,(¢/R), where E, is the secant 
modulus. 


SYMBOLS 


E, = modulus of elasticity, kilopounds per sq.in., reduced 
beyond the elastic range to compensate for plastic 
action of the material 

= secant modulus, kilopounds per sq.in. 

= panel length, in. 

= radius of curvature, in. 

= dimensionless constant 

= panel width, measured along arc, in. 

== panel thickness, in. 

= buckling strain 

fer. = buckling stress, kilopounds per sq.in. 
= ultimate stress, kilopounds per sq.in. 


INTRODUCTION 


I 1943, A SERIES OF TESTS on small curved panels in 
axial compression was conducted in the laboratories 
of The Dow Chemical Company. The results of these 
tests indicated that the effective modulus for buckling 
of such panels beyond the elastic range was closely 
approximated by the secant modulus for the panel ma- 
terial. That the secant modulus is close to the effec- 
tive modulus for buckling of flat plates has already 
been demonstrated by Gerard,! but its use in the case 
of curved plates has not, to the author’s knowledge, 
been substantiated in any published work. Since the 
completion of the 1943 tests, two additional series of 
somewhat larger panels have also been tested, and the 
results tend to bear out the earlier conclusion. This 
paper has therefore been prepared to present the data 
and show the application of the secant modulus for 
predicting buckling laads. 

For convenience, the 1943 tests are referred to as test 
series 1 and the later tests as test series 2 and 3. 


SPECIMENS AND METHOD OF TESTING 


In order to secure a wide variation in sheet proper- 
ties, panels were constructed of several different mag- 
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nesium base alloys. The longitudinal compressive 
yield stresses for the various alloys chosen are listed 
in Table 1.° Because of the variety of materials em- 
ployed, individual stress-strain curves are not given. 
The modulus of elasticity is 6,500 kilopounds per sq.in. 
in all cases. 

The specimens were rolled to curvature on a three- 
roll bending machine. In series 1, the radius of curva- 
ture was checked against a series of circular arcs drawn 
with '!/,-in. increments in radius. In series 2 and 3, 
proper curvature was obtained by using the templates 
in the loading jig.as standards. The dimensions of all 
the specimens are given in Table 2, which also presents 
the test results. The number preceding the dash in the 
specimen number indicates the test series. 

In series 2 and 3, various aspect ratios were used. 
Because the great majority of the tests had aspect 
ratios near unity, however, only such panels were in- 
cluded in the present paper. 

In the interest of simplicity, the testing methods em- 
ployed in each of the three test series will not be dis- 
cussed in detail. In general, -it can be said that the 
testing machines used were accurate within 1 per cent 
of the recorded load and that special care was exer- 
cised to secure specimens with their ends parallel and to 
obtain accurate alignment and uniform load distribu- 


TABLE 1 
Longitudinal Compressive Yield Stresses for Sheet Material 


Compressive 
Yield Stress 


Test . Thickness Kilopounds 
Series Material (In.) per Sq.In. 
1 Mh 0.025 . 24.4 

0.032 20.5 

FSa 0.016 19.4 

0.025 19.4 

0.032 19.0 

FSh 0.016 25.2 

0.025 29.3 

0.032 30.8 

J-1h 0.016 29.0 

0.025 29.2 

0.032 30.2 

2 J-ih 0.040 26.4 
0.065 23.3 

0.128 25.8 

0.248 24.6 

3 Ma 0.090 13.0 
FS-la 0.098 15.2 

FS-ih 0.098 30.2 

J-ih 0.041 31.0 

0.086 29.2 

0.092 29.2 
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1-1 
1-2 
1-3 - 
1-4 
1-5 
1-6 
ive 1-7 
ted 1-8 
>m- 1-9 
ren, 1-10 
1-11 
1-12 
1-13 
ree- 1-14 
‘Va- 1-15 
1-16 
1-17 
1-18 
ates 1-19 
all 1-20 
ents 1-21 
the 1-22 
1-23 
1-24 
sed. 1-25 
pect 1-26 
in- 1-27 
1-28 
1-29 
1-30 
dis- 1-31 
the 1-32 
cent 1-33 
1-34 
1-35 
d to 1-36 
ibu- 1-37 
1-38 
1-39 
1-40 
erial 1-41 
1-42 
ive 1-43 
ids 1-44 
n. 1-45 
1-46 
1-47 
1-48 
1-49 


Alloy 
Mh 


FSa 


FSh 


J-1ht 


R 
(In.) 
31/5 
43/, 
61/2 
63/s 
31/4 


t 


(In.) 


0 
0 
0. 
0 
0 


.0250 
.0250 
0250 
.0250 
.0255 
.0255 
0320 
0325 
0320 
.0320 
0325 
0325 
0150 
0155 
.0155 
.0155 
0250 
0250 
0250 
0250 
.0255 
.0255 
.0310 
.0305 
0310 
0310 
0310 
0305 


015 

.0165 
O15 

.0155 
0250 
.0250 
.0250 
.0255 
.0250 
.0250 
.0325 
.0320 
.0325 
.0325 
.0320 
.0325 


.026 


.026 
.035 
.016 
.016 
.026 
.026 
.026 
.026 
.035 
.016 
.016 
.026 


026 


TABLE 2 


Specimen Dimensions and Test Data 


b 
(In.) 
5.02 


(In.) 


3.20 
5.00 
5.00 


0.995 


t/R 
0.00769 
0.00769 
0.00571 
0.00571 
0.00392 
0.00392 
0.00984 
0.01000 
0.00731 
0.00711 
0.00481 
0.00473 
0.00286 
0.00295 
0.00194 
0.00221 
0.00833 
0.00833 
0.00625 
0.00625 
0.00408 
0.00408 
0.01033 
0.01017 
0.00775 
0.00775 
0.00551 
0.00530 
0.00250 
0.00275 
0.00333 
00326 
00833 
00833 
00625 
00618 
00417 
.00385 
.01083 
.01067 
.00787 
.00787 
.00533 
.00553 
0.00867 
0.00867 
0.01167 
0.01167 
0.00400 
0.00400 
0.00650 
0.00650 
0.00650 
0.00650 
0.00875 
0.00875 
0.00267 
0.00267 
0.00433 
0.00433 


~ 


pounds 
per 


18. 
22. 


bo 


(Table 2 continued on page 20) 


Cer. 
0.00346 1.072 
0.00323 1.000 
0.00272 1.134 
0.00244 1.018 
0.00194 1.179 
0.00178 1.081 
0.00453 1.097 
0.00431 1.026 
0.00355 1.156 
0.00361 1.210 
0.00237 1.173 
0.00234 1.178 
0.00106 0.891 
0.00100 0.807 
0.00091 1.1217 
0.00120 1.292 
0.00328 0.936 
0.00360 1.030 
0.00221 0.841 
0.00249 0.948 
0.00204 1.190 
0.00168 0.979 
0.00318 0.732 
0.00321 0.751 
0.00346 1.063 
0.00321 0.985 
0.00266 1.150 
0.00162 0.727 
0.00171 1.630 
0.00158 1.369 
0.00177 1.266 
0.00148 1.081 
0.00371 1.060 
0.00377 1.079 
0.00282 1.074 
0.00294 1.133 
0.00199 1.137 
0.00186 1.150 
0.00502 1.105 
0.00506 1.130 
0.00407 1.231 
0.00340 1.029 
0.00283 1.265 
0.00289 1.245 
0.00311 0.854 
0.00284 0.780 
0.00400 ~=0,816 
0.00453 0.925 
0.00155 * 0.921 
0.00170 1.012 
0.00235 0.860 
0.00253 0.925 
0:00251 0.919 
0.00289 . 1.060 
0.00368 1.001 
0.00359 0.976 
0.00118 1.053 
0.00123 1.098 
0.00182 0.999 
0.00197 1.082 


Stress 
Region * 


? 


Kilo- 
7 = a 
5.00 0.996 17.88 
4.99 5.00 1.002 17.1 
5.01 5.01 1.000 15.1 
: 5.01 5.01 1.000 13.9 
5.00 5.01 1.002 11.53 
5.02 5.01 0.998 10.7 
( 6.41 6.41 1.000 18.72 
6.43 6.40 0.995 18.3 
6.41 6.41 1.000 16.5 
0 6.43 6.40 0.995 16.65 
63/, 0 6.42 6.41 0.998 12.5 3 
67/s 0 6.42 6.41 0.998 12.4 
|| 51/4 0 3.22 3.22 1.000 6.93 
51/4 0 3.22 3.22 1.000 6.51 . 
8 0 3.23 3.22 0.997 5.93 a 
7 0 3.23 3.21 0.994 7.84 @ 
3 0 5.01 5.00 0.998 16.78 a. 
3 0 5.02 5.00 0.996 17.5 
4 0 5.02 5.01 0.998 12.9 
4 0 5.01 5.00 0.998 14.3 a 
61/, 0 6.01 5.01 1.000 12.1 
61/, 0 i 5.00 5.00 1.000 10.3 
3 0 6.42 6.40 0.997 16.9 a 
3 0 6.43 6.41 0.997 17.0 — 
4 0 6.44 6.40 0.994 17.2 = 
4 0 6.43 6.40 0.995 17.0 — 
55/s 0 6.43 6.41 0.997 15.0 _ 
53/4 0 6.43 6.41 0.997 10.0 ae 
6 0 3.21 3.21 1.000 10.1 
6 0 3.22 3.21 0.997 9.44 hs 
41/, 0 3.22 3.21 0.997 10.35 
43/, 0 3.22 3.21 0.997 8.91 . 
3 0 5.02 5.01 0.998 a 
3 0 5.00 5.00 1.000 ee 
4 0 5.01 5.00 0.998 
41/, 0 5.01 5.00 0.998 
6 0 5.01 5.00 0.998 
61/, 0 5.02 5.01 0.998 
3 0 6.42 6.41 0.998 } ; 
3 0 6.42 6.41 0.998 a 
0 6.42 6.40 0.997 
0 6.43 6.41 0.997 
6 0 6.42 6.41 0.998 
57/s 0 6.43 6.40 
| 5.00 5.00 
0 5.00 5.00 
0 6.40 6.40 
0 6.40 6.40 
0 3.20 3.20 
1-50 0 3.20 3.20 
1-51 0 5.00 5.00 : 
1-52 0 5.00 5.00 
1-53 0 5.00 5.00 = 
1-54 0 5.00 5.00 
1-55 0 6.40 6.40 : 
56 0 6.40 6.40 
1-57 0 3.20 3.20 a 
1-58 0 3.20 
1-59 0 5.00 | : 
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TABLE 2 (Continued) 
Specimen Dimensions and Test Data 


Kilo- 
pounds 
Specimen R t b L per Cer. Stress 
No. Alloy (In.) (In.) (In.) (In.) L/b t/R Sq.In. Cer. 0.42(t/R) Region* 
1-61 6 0.026 5.00 5.00 1 0.00433 11.8 0.00181 0.994 E 
1-62 6 0.026 5.00 5.00 1 0.00433 9.78 0.00150 0.824 E 
1-63 6 0.035 6.40 6.40 1 0.00583 18.2 0.00282 1.151 P, 
1-64 6 0.035 6.40 6.40 1 0.00583 18.2 0.00282 1.151 ar 
2-1 J-ih 15 0.040 14.5 18.0 1.241 0.00267 7.80 0.00120 1.070 E 
2-2 15 0.040 14.5 18.0 1.241 0.00267 7.76 0.00119 1.061 E 
2-3 15 0.065 14.5 18.0 1.241 0.00433 14.14 0.00218 1.199 E 
2-4 15 0.065 14.5 18.0 1.241 0.00433 13.27 0.00204 2.392; E 
2-5 30 0.065 6.5 6.5 1.000 0.00217 8.72 0.00134 1.470 E 
2-6 30 0.127 6.5 6.5 1.000 0.00423 15.95 0.00261 1.470 E 
2-7 15 0.128 14.5 18.0 1.241 0.00853 23.4 0.00400 1.117 P 
2-8 15 0.130 oe 1.000 0.00867 24.1 0.00418 1.149 P. 
2-9 15 0.130 14.5 18.0 1.241 0.00867 24.2 0.00420 1.153 Ff, 
2-10 60 0.248 29.8 40.2 1.349 0.00413 13.5 0.00208 1.199 E 
3-1 Ma 20.57 0.090 25.81 21.60 0.836 0.00438 10.3 0.00260 1.415 P 
3-2 20.57 0.0885 26.00 21.60 0.830 0.00430 10.4 0.00268 1.485 Pp 
3-3 15.00 0.090 25.94 25.00 0.963 0.00600 11.5 0.00360 1.430 PF 
3-4 FS-la 20.57 0.098 25.81 21.60 0.836 0.00476 11.35 0.00190 0.950 E 
3-5 20.57 0.098 26.00 21.60 0.830 0.00476 11.8 0.00205 1.026 E 
3-6 FS-1h 20.57 0.095 25.81 21.60 0.836 0.00462 14.6 0.00228 1.175 E 
3-7 20.57 0.097 25.81 21.60 0.836 0.00471 13.8 0.00218 1.102 E 
3-8 20.57 0.0985 26.00 21.60 0.830 0.00479 15.5 0.00245 1.218 if 
3-9 J-lh 20.57 0.0921 25.81 21.60 0.836 0.00448 14.07 0.00217 1.153 E 
3-10 20.57 0.0870 26.00 21.60 0.830 0.00423 14.1 0.00217 1.222 E 
3-11 20.57 0.0855 26.00 21.60 0.830 0.00416 14.9 0.00230 1.317 E 
3-12 15.00 0.0860 25.94 25.00 0.963 0.00573 14.55 0.00224 0.930 E 
3-13 15.00 0.041 25.94 25.00 0.963 0.00273 6.26 0.00096 0.837 E 


* E = within elastic range (elastic); P = beyond elastic range (plastic). 


t Only nominal dimensions recorded. 


tion in the testing machine. The methods of edge sup- 
port employed for the three series are represented in 
Fig. 1. The buckling load was determined simply by 
observation or by the occurrence of sudden “snap” 
buckles. 


RESULTS 


Test results are given in Table 2. The critical stress 
in Table 2 was determined simply as critical load/Odt. 
The critical strain was obtained from the appropriate 
stress-strain curve. The use of the values of e,+ 
[0.42(¢/R)] and thé stress region in which buckling 
occurred will be explained in the section entitled ‘“‘Dis- 
cussion.’’ Values of ultimate stress are not listed; in 
all except two instances, ultimate stress and buckling 
stress were identical. The two exceptions are as fol- 
lows: 


Specimen (Kilopounds (Kilopounds 
No. per Sq.In.) per Sq.In.) 
2-5 8.72 8.90 
2-6 15.95 16.95 


The complete test results are also plotted in Fig. 2. ° 


DISCUSSION 


A wide variety of formulas have been proposed for 
predicting the buckling stress of curved sheet in com- 
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pression. Probably the type of formula in widest use is 
for. = kE,(t/R) 


In this formula, both the constant & and the reduction 
of E, beyond the elastic range have generally been de- 
termined empirically. The tests of series 1 gave indica- 
tion that the reduced modulus could be well represented 
by the secant modulus for the panel material, or 


for, = RE,(t/R) 


Since e., = fer/E;, this latter formula can also be 
written 


Cer = R(t/R) 


Test results, 


If the secant modulus does indeed serve as a proper re- 
duced modulus, then the foregoing indicates that a 
plot of e,,, vs. t/R should give a straight line. Consider- 
ing that curved-panel data inevitably show large scat- 
ter, Fig. 2 indicates that the straight line with k = 0.42 
does provide a reasonable fit for the data. The value 
of 0.42 was chosen to make the line a lower envelope 
for the region of greatest density of the test points. It 
will be seen that the resulting straight line is in the 
lower one-third of the general scatter band. 


As evidence that the curve given by 


fer, = O.A2E,(t/R) or = 0.42(t/R) 


“4 
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is as applicable beyond the elastic range as it is within 
the elastic range, the value of e,,,/[0.42(¢/R)] was com- 
puted for each test (see Table 2). A survey of these 
values makes possible the following tabulation: 


Average Per Cent 
Value of of Points 
No. of ze Ser. on or Above 
Tests 0.42(t/R) Line 
Within elastic range 39 1.101 69 
Beyond elastic range 48 1.073 73 
All tests 87 1.085 71 


Considering the amount of general scatter, this tabu- 
lation shows conclusively that the use of the secant 
modulus for stresses beyond the elastic range predicts 
such stresses with accuracy equal to that obtained 
within the elastic range. Specifically, beyond the elas- 
tic range, relatively more points lie above the calculated 
value, but the average amount above is slightly less 
than within the elastic range. 

The value of k = 0.42 in the formula 


fers = 0.42E,(t/R) 


is not necessarily recommended for design use, since it 
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is largely up to the designer himself to decide just how 
conservative should be the formula he selects. For this 
reason, it is of interest also to note the values of k 
corresponding to the upper and lower scatter limits in 
Fig. 2. These are k = 0.60 for the upper limit and k = 
0.34 for the lower limit. This lower limit in particular 
might in some cases be selected for design use. 


CONCLUSION 


The buckling stress of curved magnesium-alloy 
panels in compression, with aspect ratios near unity 
and radius-thickness ratios between 86 and 515, can be 
predicted with equivalent accuracy within and beyond 
the elastic range by the formula 


- fer, = 0.42E,(t/R) 
where £, is the secant inodulus. 


REFERENCE 


1 Gerard, George, Secant Modulus Method for Determining Plate 
Instability Above the Proportional Limit, Journal of the Aeronau- 
tical Sciences, Vol. 13, No. 1, p. 38, January, 1946. 


Letter to 


Dear Sir: 

In connection with Conrad C. Wan’s paper, “‘Face Buckling 
and Core Strength Requirements in Sandwich Construction” 
(JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 14, No. 9, p. 
531, September, 1947), the following comments might be of in- 
terest. It is known that no structural element ever fails because 
of instability but only because the stresses in some part of it at- 
tain too high values. Instability phenomena can bring about 
such excessive stresses, but, in the end, failure is due to the ex- 
cessive stress and not directly to instability. For this reason it is 
gratifying to see that the problem of the instability failure of 
sandwich panels was approached by the author from the stress 
standpoint. 

This approach leads to the inevitable conclusion that the 
buckling stress of the faces depends upon the ultimate tensile or 
compressive stress of the core material, as well as upon its elastic 
modulus. In contrast, pure stability calculations yield formulas 
from which the ultimate stress is missing. Indeed, it is one of the 
most important conclusions drawn by Mr. Wan in his paper that 
the ultimate stress of the core is a decisive factor in face buckling. 
Yet his own experiments contradict him, as may be seen from 
Fig. 5 of the paper. The experimental failing stress values of the 
sandwich face vary between 64,000 and 76,000 Ibs. per sq.in., 
and the averages plot almost exactly along a horizontal line cor- 
responding to 70,000 Ibs. per sq.in. This means that, in the range 
of core strength investigated—that is, when the ultimate tensile 
stress of the core varies from 750 to 1,750 lbs. per sq.in.—the fail- 
ing stress of the face is independent of the ultimate stress of the 
core. 

It may be mentioned that simple columns fail slightly below the 
classical instability limit when manufactured and centered with 


the Editor 


reasonable care, while the failing load of thin plates supported 
along all the four edges is much higher than their stability limit. 
Thin-walled circular cylinders fail well below the classical buc- 
kling load. The absence of a core strength influence upon the face 
failing load in Mr. Wan’s tests might indicate that in this case 
again the classical buckling load is not too different from the ac- 
tual failing load. 

The writer is perfectly willing to admit that the cube root 
formula derived by him jointly with Dr. S. E. Mautner is just 
an approximation. Nevertheless, it does not seem to fare so 
badly in the example cited by Mr. Wan. A reduction of the 
buckling stress from 212,000 to 70,000 Ibs. per sq.in. would be 
possible only if the reduced modulus of 75S-T Alclad were as- 
sumed as 400,000 Ibs. per sq.in., which is an extremely low value. 
On the other hand, no reliable modulus values can be had for 
strains of such magnitude. When the buckling stress calculated 
with a reduced modulus is higher than the yield point stress, the 
writer usually assumes the yield point stress as the maximum that 
can be supported by the material. According to the Aluminum 
Company of America, the yield point stress of 75S-T Alclad is 
66,000 Ibs. per sq.in. and its ultimate stress is 77,000 lbs. per 
sq.in. It is of interest to note that all specimens shown in Fig. 5 
failed between 64,000 and 76,000 Ibs. per sq.in. 

Finally, the writer wishes to express his hope that Mr. Wan 
will continue his valuable analysis in order to clarify the reasons 
for the discrepancy between theory and experiment. 


N. J. Horr 
Professor of Aeronautical Engineering 
Polytechnic Institute of Brooklyn 
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Prediction and Reduction to Minimum 
Properties of Plate Compressive Curves 


E. H. SCHUETTE* anp J. C. McDONALD* 
The Dow Chemical Company 


SUMMARY 


A method is presented for predicting plate compressive curves 
or reducing them to minimum properties, based on the similari- 
ties that exist between such curves and the stress-strain curve for 
the material. Comparison with experimental results indicates 
that the method is satisfactory for extruded material. There are 
no check data for sheet material; it is possible that predicted 
curves for sheet material may be unconservative, though not 
greatly so. 


INTRODUCTION 


6 le SUBJECT OF FLAT PLATE INSTABILITY in com- 
pression has received extensive experimental 
treatment by the N.A.C.A. for aluminum alloys'~ and 
has also been studied experimentally by the same group 
for one magnesium alloy.’ The results of these studies 
have been embodied in plate compressive curves for 
the various materials with more-or-less typical proper- 
ties. As it is frequently required that designs be made 
on the basis of material with minimum guaranteed 
properties, it seems desirable to have available a method 
for estimating plate compressive curves for minimum- 
property material. Moreover, an estimate of such 
curves for new or untested materials without extensive 
experimental work would evidently be of considerable 
value. 

It has already been demonstrated by Gerard® that 
the effective modulus for buckling of flat plates in 
compression is closely approximated by the secant 
modulus. That is, the shape of the curve of experimen- 
tal buckling stress against calculated elastic buckling 
strain approximates the shape of the stress-strain curve. 
In the present paper, it is pointed out that a rather 
consistent relationship between the buckling curve and 
the stress-strain curve exists for extruded material, dif- 
ferent from that shown by, sheet material. It is further ~ 
shown that the maximum average stress is related to 
the stress-strain curve in a fairly consistent fashion. 
On the basis of such relationships, a method is outlined 
for predicting plate compressive curves for any material 
with a known stress-strain curve, and an example is 
given to show the accuracy obtained. 


SYMBOLS 
E = compressive modulus of elasticity, kilopounds per 
sq.in. 


Received May 19, 1947. , 
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b/t = width-thickness ratio of plate 
bw/t = width-thickness ratio of web (references 1-7) 


k = dimensionless coefficient in the formula for plate 
buckling 

n = dimensionless coefficient for reduction in E beyond 
the elastic range 

m = Poisson’s ratio 

o = compressive stress, kilopounds per sq.in. 

Cc.y, = compressive yield stress 

Gcr. = critical compressive stress 

Tmar. = average compressive stress at maximum load 


RELATIONSHIP BETWEEN PLATE CURVES 
AND STRESS-STRAIN CURVES 


The method for determining the buckling stress of a 
flat plate in compression when plate curves are avail- 
able, as set forth by the N.A.C.A.,' is based on the 
buckling formula 


Ger, = E/(12(1 — 
This formula is generally written in the form 
Ger,/n = kw®E/(12(1 — w*)(b/t)?) 


from which o,,,/n can be computed using the elastic 
modulus E. The value of o,,, is then determined by use 
of a curve of VS. 

Because the supported edges of the plate are pre- 
vented from undergoing lateral displacements, they re- 
main straight when the plate buckles, and the plate is 
thus enabled to carry some additional load after buck- 
ling takes place. There is therefore a maximum load 
for a plate which is higher than its critical, or buckling, 
load. The value of ¢.¢;, the average stress at maximum 
load, is generally determined by a calculation for o,,./n 
and a curve Of Gmaz, VS. Ger,/N- 

If o.;,/n is divided by £, the resulting quantity is in 
strain units; thus, curves of o,. vs. (¢¢,,/n)/E and 
maz, VS. (Gcr,/n)/E can be plotted on the same coordi- 
nate system with a stress-strain curve. When such plots 
are made on the basis of presently available test data, 
it is immediately apparent that the type of relationship 
between plate curves and stress-strain curves is different 
for extruded material from that for sheet material. A 
possible cause of the difference will subsequently be 
discussed. The relationships are typified by the curves 
in Fig. 1, prepared from data in reference 4, and Fig. 2, 
prepared from data in reference 2. The additional in- 
formation on Fig. 2 will be discussed in the section that 
follows. 
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Fic. 1. Comparison of plate compressive curves and stress- 
strain curve for extruded 24S-T aluminum alloy (¢,.y, = 49.9 kilo- 
pounds per sq.in.), typical of extruded material. 
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Fic. 2. Comparison of plate compressive curves and stress- 
strain curve for 17S-T aluminum-alloy sheet (0.4. = 41.4 kilo- 
pounds per sq.in.), typical of sheet material. | 


TABLE 1 


Relationships Between Plate Compressive Curves and Stress-Strain Curve for Aluminum-Alloy Sheet 
(o¢.y. = 41.4 kilopounds per sq.in. for 17S-T, 45.1 kilopounds per sq.in. for 24S-T Z’s, and 43.6 kilopounds per sq.in. for 24S-T channels) 


o1/00— o2/a0 
Offset from 24S-T 2 24S-T 24S-T 24S-T 
Modulus Line 17S-T (2) (channel) Average 17S-T (Z) (channel) Average 
0.001 1.048 1.024 1.060 1.044 1.130 1.100 1.132 1.121 
0.002 1.090 1.077 1.115 1.094 1.150 1.142 1.175 1.156 
0.003 1.086 1.088 1.112 1.095 1.132 1.137 1.158 1.142 
- €2/€0 
Fraction of 24S-T 24S-T 24S-T 24S-T 
Yield Stress 17S-T (2) (channel) Average 17S-T (Z) (channel) Average 
0.45 1.000 1.000 1.000 dike 0.317 0.292 0.264T 
0.50 1.015 1.000 1.008 | 0.408 0.402 0.351T 
0.55 pas 1.022 1.017 1.020 es 0.520 0.542 0.460T 
0.60 1.035 1.029 1.032 0.685 0.749 0.621f 
0.65 1.022 1.040 1.034 1.032 0.562 0.938 0.947 0 816 
0.70 1.035 1.050 1.041 1.042 0.726 0.988 0.989 0.991 
0.75 "1.043 0.879 0.985 * 


* No averages taken. 


+ Estimated on the assumption that average of all three would have same relation to average of last two only as it has for 0.65 o¢.y. 


PREDICTION OF PLATE CURVES FROM 
STRESS-STRAIN CURVES 


The specific relationships used in predicting the plate 
curves are illustrated by Fig. 2. The stress-strain curve 
is an average of curves B and C in Fig. 4 of reference 2 
(Table I of reference 2 shows that all with-grain Z and 
channel columns were cut from the sheets represented 
by curves B and C). The curve of o¢,. vs. (¢¢,/n)/E 
was obtained by dividing the abscissa values in Fig. 9 
of reference 2 by the modulus of elasticity, 10,700 
kilopounds per sq.in., and the curves of @maz, VS. 


(oer./n)/E was obtained by performing the same divi- 
sion on abscissa values in Fig. 11 of reference 2. In the 
latter case, the left end of the curve is that for b,/t = 
32, and an average of the curves for Z’s and channels 
was used. The highest value of },,/t for which there 
were test data was selected in order to minimize the ef- 
fect of the forming operation on material properties (see 
discussion of this point in reference 1, page §). 
Referring to Fig. 2, let us designate values of stress 


and strain on the stress-strain curve by oo and 6, re- 
spectively. Similarly, values of ordinate and abscissa 
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TABLE 2 


Relationships Between Plate Compressive Curves and Stress-Strain Curve for Extruded Aluminum Alloys 
* (cey-. = 49.9 kilopounds per sq.in. for 24S-T, 58.6 kilopounds per sq-in. - 72.1 kilopounds per sq.in. for R303-T, and 79.2 kilopounds per sq.in. for 


Offset 
from 
Modulus 
Line 24S-T 14S-T R303-T 75S-T Average 24S-T 14S-T R303-T 75S-T Average 
0.001 0.923 0.895 0.941 0.914 0.918 0.979 0.948 0.980 0.959 0.966 
0.002 0.944 0.922 0.957 0.927 0.938 0.988 0.967 0.993 0.963 0.978 
0.003 0.955 0.938 0.968 0.932 0.948 0.992 0.983 1.001 0.964 0.985 
0.004 0.965 0.955 0.975 0.931 0.956 1.010 0.999 1.008 0.962 0.995 
Fraction 
of Yield = «/e 
Stress 24S-T 14S-T R303-T 75S-T Average 24S-T 14S-T R303-T 75S-T Average 
0.60 ose 1.000 1.003 1.002 0.437 0.529 0.416* 
0.65 1.022 1.008 1.004 1.008 1.010 0.552 0.443 0.642 0.738 0.594 
0.70 1.024 1.030 1.017 1.014 1.021 0.720 0.735 0.831 0.909 0.799 
0.75 1.039 1.058 1.035 1.023 1.039 0.955 0.913 0.947 0.982 0.949 
0.80 1.055 1.101 1.049 1.045 1.062 0.991 0.998 1.009 1,008 1.002 


* Estimated on the assumption that average of all four would have same relation to average of last two only as it has for 0.65¢¢.y.. 


on the curve of oe, vs. (¢¢r,/n)/E are designated o; and 
€, and values on the curve of &maz. VS. (Gcr./n)/E are 
designated o2 and €. Also included on the plot (Fig. 2) 
are sloping lines representing various offsets from the 
modulus line and horizontal lines representing various 
decimal parts of the yield stress. For a given offset, one 
can pick off values of o0, 01, and o2 and compute the 
ratios o;/o) and o2/o. Also, for a given decimal part 
of the yield stress, one can pick off €, 4, and € and com- 
pute the ratios €,/¢ and 

Curves similar to those in Figs. 1 and 2 were also 
plotted for 24S-T sheet,’ extruded 14S-T,*® extruded 
R303-T,° and extruded 75S-T,* and a similar procedure 
was followed in each case. The results are tabulated 
in Tables 1 and 2. For those extruded materials where 
Z’s and channels produced curves slightly different 
from those obtained with H-sections, an average curve 
was used. Presumably, as is pointed out in the refer- 
ences, this difference is largely due to the fact that 
machining off flanges from an extrusion of somewhat 
variable properties produces a section whose average 
compressive properties are not the same as those of the 
original section. On the basis of this explanation, the 
use of an average curve seems desirable. 


The curves for extruded O-1LHTA magnesium alloy’ 
were not available at the time the original work was 
done. Because this original work was completed with- 
out the O-1HTA data and because the present existence 
of those data furnishes an excellent means for checking 
the accuracy of the method, they were not included in 
the preparation of Tables 1 and 2. 

A brief study of Tables 1 and 2 will serve to show that 
two distinct relationships do exist between plate curves 
and stress-strain curves, one for extruded material and 
one for sheet material; the greatest distinction is evi- 
dent in the values of the stress ratios for the various 
offsets from the modulus line. 


50 


40 


002-4 € 


Fic. 3. Stress-strain curves used in predicting plate curves for 
magnesium alloys of minimum guaranteed properties. 


Averages of the several ratios are listed in Tables 1 
and 2. To predict the plate curves for any material 
whose stress-strain curve is given, it is only necessary 
to apply the average ratio values from Table 1 or 2, 
depending on whether sheet or extruded material is 
under consideration, to produce points on the desired 
curves of vs. and VS. (Ger,/n)/E. 
Subsequent multiplication of abscissa values by E 
gives the curves in their final form. 

As an illustration of the method, plate curves were 
prepared for several magnesium alloys (FS-1h,O-1HTA, 
and two experimental materials, ZK60 and ZK60A) 
with minimum guaranteed properties. The stress- 
strain curves on which the work was based are shown 
in Fig. 3, and the resulting plate compressive curves are 
presented in Figs. 4-7. 


ACCURACY OF METHOD 


In judging the accuracy of the method just presented, 
it must be remembered that all the test data on which 
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Fic. 4. Predicted plate curves for FS-1h magnesium-alloy 
sheet of minimum guaranteed properties (o.y, = 25 kilopounds 


per sq.in.). 
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Fic. 5. . Predicted plate curves for extruded O-1HTA magne- 
sium alloy of minimum guaranteed properties (o-.y. = 30 kilo- 
pounds per sq.in.). 
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Fic.6. Predicted plate curves for extruded ZK60 experimental 
magnesium alloy of minimum guaranteed properties (oy. = 
30 kilopounds per sq.in.). 
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Fic. 7. Predicted plate curves for extruded ZK60A experi- 
mental magnesium alloy of minimum guaranteed properties 
(¢y.c. = 35 kilopounds per sq.in.). 
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PREDICTION OF PLATE COMPRESSIVE CURVES 


the curves are based were obtained from H-, Z-, or 
channel sections. There is little known that would 
enable one to predict how well such curves will describe 
the action of other kinds of plate assemblies. The 
method is intended simply as a means of arriving, with 


few or no tests, at the same curves that would be ob-. 


tained from tests of H-, Z-, or channel sections.. The 
only applicable measure of its accuracy, therefore, is 
how well it will predict just that type of curve. 

As previously stated, reference 7, presenting experi- 
mentally determined curves for O-1HTA, was pub- 
lished subsequent to the original work on predicting 
curves. As a check on the method, curves were then 


predicted on the basis of the stress-strain curve for the . 


material used in reference 7 (average of two curves for 
local-instability columns, Fig. 5, reference 7), and 
these curves wére compared with the N.A.C.A. test 
curves. The comparison is shown in Fig. 8; it pro- 
vides satisfactory substantiation of the procedure. The 
maximum unconservatism evidenced by the predicted 
curves is of the order of 6 to 7 per cent, in the region of 
Ser, 15 kilopounds per sq.in. In this same region, 
the overall test scatter band’ is of the order of 20 per 
cent. 


DISCUSSION 


Predicted curves for extruded material have been 
satisfactorily checked by test experience, as evidenced 
by Fig. 8. There is admittedly still room for doubt in 
the case of sheet material. The first question that 
arises in this connection is why the curves for sheet ma- 
terial, in the high-stress region, consistently run so 
much higher relative to the stress-strain curve than do 
those for extruded material. It seems unlikely that 
this difference is due to the greater stability of the 
curved corners, compared to square corners in extru- 
sions, when the material is tested in the form of bent-up 
Z’s and channels; for, if the curved corners did appre- 
ciably increase stability, the effect should also be 
noticeable in buckling stresses within the elastic range. 
In references 1 and 2, no such effect is apparent in the 
low-stress region. It is known, however, that the form- 
ing operation does appreciably increase the compres- 
sive yield stress in the corners. Possibly, then, a higher 
stress-strain curve should be used in conjunction with 
plate curves based on tests of bent-up sections for ma- 
terials where bending causes an increase in yield stress, 
but it is difficult to prescribe the proper stress-strain 
curve. 

If the forming operation is indeed responsible for* 
raising the plate curves, then the curves predicted for 
FS-1h sheet are probably slightly unconservative in the 
high-stress region because no increase in corner proper- 
ties takes place in the hot-forming of FS-lh. More- 
over, the curves of references 1 and 2 will be unconserv- 
ative if applied to plates where no forming is involved. 
Further valuable information on this subject could be 
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Fic. 8. Comparison of predicted and experimental (N.A.C.A.) 
plate curves for extruded O-1HTA magnesium alloy (oy. = 
33.6 kilopounds per sq.in.). 


obtained from tests of sheet material with uniform prop- 
erties in the bent-up condition, such as FS-1h sheet. 
The method has been presented from a standpoint 
of predicting curves for materials for which yo test 
data are available. If it is desired to reduce curves al- 
ready obtained experimentally to minimum guaranteed 
properties, it would probably be best not to use the 
average values of the stress and strain ratios given in 
Tables 1 and 2 but rather to use the specific values ob- 
tained from the test curves for the material considered 
and apply these to a minimum stress-strain curve. 


CONCLUSION 


Comparison with experimental results indicates that 
the method presented for predicting plate compressive 
curves or reducing them to minimum properties is satis- 
factory for extruded material. There are no check 
data for sheet material; it is possible that predicted 
curves for sheet material may be unconservative, 
though not greatly so. 
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A Method of Wind-Tunnel Testing Through 


the Transonic Range 


JOHN H: WEAVER* 
Lockheed Aircraft Corporation 


SUMMARY 


Wind-tunnel testing by conventional methods in the so-called 
transonic range is extremely difficult and unreliable because of 
the blocking phenomenon between the model and the tunnel 
walls. This paper outlines a method of modifying existing high- 
speed wind-tunnel test sections with a suitably contoured bump 
on which reflection plane models can be satisfactorily tested 
through the speed of sound. The method depends upon the es- 
tablishment of a local high-velocity region over the bump in 
which half-models mounted vertically can be tested. Experi- 
mental data on several models tested by this method are pre- 
sented, together with typical velocity distributions over the 


bump contour. 
The results of the tests made to date appear to be reliable, and 


the method offers a means of obtaining rapidly, and at moderate 
costs, large amounts of aerodynamic data in the all-important 
region near M = 1.0. 

INTRODUCTION 


| rms OF THE CHOKING PHENOMENON between the 

model and the air-stream boundaries, conventional 
wind tunnels cannot be used to obtain continuous data 
through the speed of sound. Testing can be done be- 
low a Mach Number of 1 in a subsonic tunnel or above 
a Mach Number of 1 in a supersonic tunnel, leaving a 
“forbidden” band between the two the width of which 
is determined by the care and diligence exercised in 
conducting the tests. Unfortunately, large changes in 
aerodynamic forces are likely to occur on the airplane 
at or near M = 1.0, and, consequently, extrapolation of 
data through this range is difficult and unreliable. It 
is imperative that a large amount of experimental data 
be obtained near M = 1.0 before airplanes can be 
designed to fly safely at these speeds. Some success 
has been realized by testing with one of the following 
two methods: (1) free air flights of models or missiles! 
and (2) flight tests of reflection plane models mounted 
on the wings of fast-flying aircraft.2 Both of these 
methods are costly and not suitable for extensive re- 
search work, since they will not yield large amounts of 
data on numerous configurations in short periods of 
time. 


The technique reported here consists of installing a ° 


suitably contoured bump in the test section of a con- 
ventional subsonic wind tunnel so that a supersonic 
velocity region is obtained over the bump in which 


Presented at the Annual Summer Meeting, I.A.S., Los Angeles, 
August 7-8, 1947. 
*Department Manager, Wind Tunnel. Now with Telecom- 


puting Corporation. 


half-span models mounted vertically can be tested. 
This work was started in August, 1946, and completed 
in June, 1947, by Lockheed Aerodynamic Laboratory 
personnel. The bump was designed for use in the 
Southern California Cooperative Wind Tunnel (CWT) 
in Pasadena, Calif. 


APPARATUS AND METHODS 
Bump Configuration 


A diagrammatic sketch of the bump is shown in Fig. 
1, and a photo of the installation in CWT is given in 
Fig. 2. In order to develop a suitable configuration 
with a minimum expenditure of time and money, a 
1/9 scale model of CWT was constructed consisting of 


TEST SECTION CaRT ——— 


\ 
\——REFERENCE STATC ORIFICE 


TOTAL HEAD TUBES 


fy —BOWNDARY LAYER DUCT 
p 


Fic. 1. Sketch of bump installation in CWT. 


Fic. 2. Photo of bump installation in CWT. 
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WIND-TUNNEL T,ESTING THROUGH THE TRANSONIC RANGE 


Fic. 3. Photo of !/22 scale model of CWT. 


an entrance cone, test section, diffuser, and centrifugal 
blower unit, as shown in Fig. 3. The bump contour 
was developed in this model tunnel by cut-and-try 
processes, with the following practical factors greatly 
influencing its overall size and shape: 


(1) The bump should be no longer than 15 ft. so as 
to be entirely confined to the removable cart forming 
the tunnel test section. 


(2) The bump should be thick enough to completely 
house a strain-gage balance system and remotely con- 
trolled angle-changing apparatus. 

(3) The models to be tested would have 5- to 8-in. 
semispan wings. 


Preliminary tests in the model tunnel indicated that 
the magnitude of the local Mach Number and the ver- 
tical velocity gradient above the bump contour were 
controlled primarily by the ratio of bump thickness to 
tunnel height. Obtaining small velocity gradients 
both vertically and horizontally was considered more 
important than high testing speeds; consequently, 
the final bump contour was obtained by modifying a 
circular are contour of proper thickness to give a 
maximum local Mach Number of about 1.2 at its center 
line. The thick boundary layer next to the tunnel 
walls was ducted under the bump so as not to interfere 
with the model tests, Ordinates for the final bump are 
listed in Table 1. 


Rather large pressure differentials were induced at 
various points in the tunnel test section because of the 
bump installation, and, consequently, the small tunnel 
tests were helpful in determining where the test section 
structure had to be strengthened to carry the bump 
loads. 


TEST WING 

COVER PLATE 

CLEARANCE 

WING MOUNT PLATE 

TURNTABLE 
BALANCE PLATE 


STRUTS 
~——DRAG STRUT 
LOAD STRUT 
BOX 
FLEXURE 
“STRUT INSULATOR (TYP) 

~ TURNTABLE. BOX 

-DRIVING MECHANISM 


Fra. 4. Sketch of six-component balance system. 
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TABLE 1 
Bump Ordinates 


Distance Distance 
Aft of Leading Above Tunnel 
Edge, In. Floor, In. 
0 2.552 
8.756 6.820 

19.756 11.220 
30.756 14.564 
39.556 16.588 
48.356 17.908 
52.756 18.348 
57.156 18.656 
61.556 18.854 
65.956 18.942 
70..356 18.942 
74.756 18.876 
79.156 18.700 
83.556 18.414 
87.956 17.952 
96.756 16.544 

105. 556 14.520 : 

Fic. 5. Flow calibration plate 

132.020 4.510 


TABLE 2 
Balance Characteristics 


Maximum 
Load, Accuracy, 
Component Lbs. Lb. 

Front lift 100 +().2 
Rear lift 50 +().1 
Drag 10 +0 .02 
Side force No. 1 40 +0.08 
Side force No. 2 40 +0.08 
Side force No. 3 80 +(0).16 


Balance System 


A six-component strain-gage balance system was 
designed to be mounted within the bump contour. A 
sketch of the installation is shown in Fig. 4. In order 
to make the system as simple and reliable as possible, 
the whole balance unit and model were rotated about a 
vertical axis for angles-of-attack changes. This meant 
that the balance measured airplane axis forces and 
moments that had to be transferred to conventional 
wind axis in the data reduction. A 16-in. diameter 
turntable was provided through which the models 
passed with small clearances so that only air loads 
on the model were measured by the balance sys- 
tem. 

Considerable trouble was encountered with the first 


balance system because of severe model vibration above 


critical speeds. A redesign with a reduction in sensi- 
tivity and better mounting of the strain gages on the 
bending beams resulted in a satisfactory system, which, 
to date, has operated for 200 hours of testing with no 
trouble. The characteristics of the balance when read 
by a Browne self-balancing potentiometer are listed in 
Table 2. 


Bump Velocity Survey 


A vertical pressure plate was installed in the model 
test region having flush orifices covering an area 10 
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Fic. 6. Mach Number distribution above bump surface. 


in, above the bump and 8 in. fore and aft of the maxi- 
mum ordinate (see Fig. 5). A flush reference static 
orifice was installed approximately 3 ft. to the side of 
the bump centerline, and two total head tubes were 
mounted forward in the subsonic flow region (see Fig. 
1). From these readings a calibration curve of average 
Mach Number at the model vs. Mach Number indi- 
cated by the reference tube was obtained. A typical 
Mach Number distribution above the bump is shown 
in Fig. 6. The boundary-layer profile at the model 
location was measured with a total head rake and is 
shown in Fig. 7. In order to account for this 0.6-in. 
thick boundary layer, the zero stations of the models 
were moved out 0.2 in. from the bump contour, and all 
coefficients were based on the model dimensions out- 
side this 0.2-in. station. 


Models 


The semispan models were constructed of solid brass 
or aluminum alloy, the former being preferred where 
possible, since the component parts could then be 
readily soldered. It was found that the model con- 


30 JOURNAL OF THE AERONAUTICAL SCIENCES—JANUARY, 1948 
FI 
Bump Bounoary ff / 71200 
LAYER 
mo 
i Tes 
dec 
was 
Ibs. 
van 
( 
plo 
test 
3 
bers 
Nut 
i 
win; 
. 


WIND-TUNNEL TESTING THROUGH THE TRANSONIC RANGE 31 


10 


| VELOCITY IN BOUNDARY LAYER 


VoxVELOCITY OUTSIDE BOUNDARY LAYER 
9-f V2 


3 
DISTANCE ABOVE BUMP SURIACE— 


Fic. 7. Boundary-layer profile at centerline bump ~ M = 
1.000. 


Fic. 8. Close-up view of model installed for test. 


tours could readily be held to +0.002-in. tolerance for 
moderate model costs. A photo of the P-80 model 
mounted for test is shown in Fig. 8. 


Testing Technique 


After some preliminary tests had been made, it was 
decided that the best operating technique to be followed 
was that of setting a constant value of '/.9V? = 500 
lbs. per sq.ft. for all Mach Numbers tested. This 
could be done conveniently at CWT because of the 
variable pressure feature of the tunnel and was ad- 
vantageous for the following reasons: 

(1) Data reduction was greatly simplified, allowing 
plots of final corrected coefficients to be made while the 
test was in progress. 

(2) The sensitivity of the balance system gave 
the same accuracy to the coefficients for all Mach Num- 
bers. 

The curve in Fig. 9 shows the variation of Reynolds 
Number vs. Mach Number using the above operating 
technique based on an average chord length for the 
wings tested. 


Test DATA 
P-80 Wing 


Figs. 10 and 11 show the effect of Mach Number on 
P-80 wing alone data at constant values of C,;. The 
large change in zero lift angle and accompanying Cy 
changes are apparently due to the airfoil section cam- 
ber causing unsymmetrical shock formation on the up- 
per and lower surfaces. Above M = 0.975, the coef- 
ficients have stabilized again with the aerodynamic 
center of the wing shifted aft to 40 per cent M.A.C. 
and the zero lift angle shifted to 0°. 


Q 6 STAGNATION TEMP + 90°F 
x AVERAGE CHORD 3 1M 
> 500 <88/cr* FOR All MACH W/O 
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Fic. 9. Average Reynolds Number vs. Mach Number during 
test. 
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Fic. 10. P-80 wing alone. 
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Fic. 11. P-80 wing alone. 
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Fic. 12. Effect of increasing and decreasing Mach Number on 
P alone. Cpa =0° 


Tests on models mounted by conventional methods 
in high-speed wind tunnels have shown considerable 
hysteresis in the coefficients as the Mach Number was 
increased or decreased above the wing critical. A simi- 
lar test was conducted on the P-80 wing on the bump 
with no effect, as shown in Fig. 12, indicating that the 
hysteresis is a tunnel-blocking phenomenon and not 
likely to occur in flight. 
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Fic. 13. Effect of leading-edge sweepback angle on Cp vs. M. 


Sweepback Effect 


Although no direct attempt was made to test a series 
of wings with sweep angle as the only variable, the data 
shown in Fig. 13 form a fairly true family of curves for 
the effect of leading-edge sweep angle. From these 
data, it is apparent that leading edge angle has four 
effects on Cp vs. M, as follows: 


(1) The knee of the drag curves move toward higher 
values of M as the leading-edge sweepback angle in- 
creases. Although no definite law can be derived from 
these data, the variation is fairly close to proportional 
to the square root of the cosine of the leading-edge 
angle. 

(2) The slope of the drag coefficient rise before the 
knee of the curve decreases with increasing leading-edge 
angle. 

(3) The slope of the drag coefficient rise after the 
knee of the curve decreases with increasing leading-edge 
angle. 

(4) The peak drag coefficient decreases with in- 
creasing leading-edge angle. 

A pronounced improvement was also noticed on the 
smoothness of the balance readings above the criti- 


cal Mach Number for wings with more leading-edge. 


angle. 

The angle of sweep of the maximum thickness line of 
wings is a factor effecting their high-speed performance, 
as shown in Fig. 14. A gain in critical Mach Number 
for all values of C, tested was realized for a given 
plan-form wing by sweeping the maximum thickness 
line. 
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Ww [Wa 
ROOT THICKNESS MPC HAC 
TIP THICKNESS Bec 


ROOT MAX. ORDINATE |40%C 
TIP MAX ORDINATE \40%C! 


ASPECT RATIO 37 |37 
TAPER RATIO 2 
| | | | | 
M 
Fic. 14. Effect of maximum thickness sweepback angle on Cp 
vs. M. 
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FUSELAGE DATA 
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CROSS SECTION SHAPE -ELIDSE 
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. 
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TAIL DATA 


ASPECT RATIO - 3.2 
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SECTION - 654008 
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Fuselage and horizontal tail drag data tested alone 
and with other model components. 


Fic. 15. 


Miscellaneous Data 


Some testing was done with fuselage and horizontal 
tails on the models. The results indicate that this 
method of testing will be useful in determining inter- 
ference effects, as shown in Fig. 15, and tail efficiency 
variations with Mach Number. Elevator effectiveness 
(or the lack of it) can be determined by building a series 
of tails with preset elevator angles. 


P.80 WING | SWEPT-BACK WING 
ASPECT RATIO 6.4 3.0 
LEADING EDGE ANGLE| /Q°- 30’ 49°-/5' 
ROOT AIRFOIL 66, | 65, -0068 
TIP AIRFOIL 65, 213, A-.5 | 65, -008 
TAPER RATIO 26 2 
GEOM TWIST o° 
= 
SWEPT-BACK WING 
6 8 3 Lo 
Fic. 16. .Spanwise center of lift variation with Mach Number 
Cy =0.10. 


Spanwise center of lift as determined from rolling 
moment of the half-span wings is shown in Fig. 16 
for the P-80 wing and a sweptback wing. This type 
of data is useful in structural design and in determining 
the spanwise region of any loss in lift with increasing 
Mach Number. 


CONCLUSIONS 


The bump testing technique offers for the first time 
a method of wind-tunnel testing through the speed of 
sound. Existing high-speed wind tunnels throughout 
the country can be suitably modified for very little cost 
so that transonic testing can be undertaken without 
restricting the normal use-of the tunnel. This will 
mean a great saving in time and cost of developing fu- 
ture military aircraft and will result in greater safety 
for the pilots who fly them. 

The question of reliability of the data cannot be 
answered quantitatively, since there are no reliable 
data in the transonic range against which to check. 
The following are some of the factors that indicate that 
the bump data will come close to predicting true flight 
characteristics: 


(1) Tests-on similar wings by conventional wind- 
tunnel methods check the bump results up through the 
critical Mach Numbers of the wings. 

(2) Noeffect of models mounted on the bump could 
be found on the walls or ceiling pressure distribution. 
As an indication of the rapidity with which the effect 
of the model dies out with distance, the data in Fig. 17 
are presented. 

(3) Data could be repeated on the same models 
after removing and reinstalling the bump several weeks 
later within +2 per cent. 
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Improvements in the method can and will be made, 
such as the following: 

(1) Reduction in boundary-layer thickness by more 
adequate removal just ahead of the models. 

(2) Obtaining higher test Mach Numbers with no 
increase in velocity gradient by the use of auxiliary 
ceiling bumps. Some work, with favorable results, has 
already been done on this by the author. 

(3) Improved shape of the aft portion of the bump 
to cut down tunnel power consumption so that higher 
tunnel pressures can be used with resulting higher test 
Reynolds Numbers. 

(4) Downstream choking of the main air stream to 
give steadier testing conditions over the bump. This 
method has proved successful for conventional tunnel 
testing.® 
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Letter to the Editor 


Dear Sir: 

In my article ‘‘Some Ballistic Contributions to Aerodynamics” 
(JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 14, No. 3, 
pp. 155-166, March, 1947), I state in a footnote on page 163: 
“Tt has been shown that the base pressure does not change with 
length for an angle of 7° 15’, and it seems reasonable to believe 
that this relation will hold at other small angles. The base 
pressure is therefore equal to the pressure on a square-based 
body and is invariant with both length and angle.” 

In view of experimental data that I have obtained subsequent 
to the writing of my article, it is my opinion that the above 


statement is too sweeping, although at the time it appeared to 
me to be a reasonable interpretation of the data obtained. Mis- 
sile designers should not assume it to be a demonstrated fact 
that the base pressure is independent of the angle and length of 
a conical boattail. * 


A. C. CHARTERS 
Ballistic Research Laboratories 
Department of the Army 
Ordnance Department 
Aberdeen Proving Ground 


pin REF STATIC 
ORIFICE 
M 
; ” An 
methc 
mechs 
4 
4 
P* 
meth 
a pre 
(b) 
(c) t 
dem 
cout 
velo 
has 
ball 
‘ is 1 
mal 
cou 
refe 
1 
rise 
frec 
bee 
she 
osc 
the 
f 
i I 
: 4 
Me Dir 
4 
‘ 
tel] 


listribu- 


"hrough 
od, Na- 


Shock- 
autical 


‘ed to 
Mis- 
| fact 
th of 


RS 

ories 
irmy 
nent 
und 


Loss of Spin of Projectiles 
Part I—Experimental Method* 


J. A. VAN ALLEN? 
Applied Physics Laboratory, The Johns Hopkins University 


ABSTRACT 


A novel method is described for the measurement of the spin 
of projectiles throughout flight. The essential element of the 
method is a small rugged radio transmitter in the nose of the 
projectile with radiation pattern not symmetrical about the 
mechanical axis. Sample data are given, and applications are 
suggested. 


INTRODUCTION 


ye TO THE WORK DESCRIBED HEREIN, the spin 
of shells in flight has been measured by: (a) the 
method of R. H. Kent,* in which a shéll, equipped with 
a projecting metal pin, is fired through a series of cards; 
(b) the photographic method of von Neesen,** and 
(c) the mechanical time fuse methed of Hill.** 

Measurements made by the above methods clearly 
demonstrated the loss of spin of shells in flight 
and yielded determinations of the spin-destroying 
couple. 

The radiosonde method of the present report was de- 
veloped concurrently with the radio-proximity fuse and 
has been shown to possess considerable power as a 
ballistic and aerodynamic technique. 


BASIS OF THE METHOD 


A small battery-powered radio-frequency oscillator 
is located in the nose of the shell in the position nor- 
mally occupied by the nose fuse. This oscillator is 
coupled to an antenna of such a character that the radia- 
tion pattern of the system is not axially symmetrical. 
The combination of oscillator, battery, and antenna is 
referred to as a radio spin sonde. 

Thus a sonde-equipped spinning shell in flight gives 
rise to an amplitude and polarization modulated radio- 
frequency signal at any fixed position. 

In practice, the radio frequency of the sonde has 
been selected so that the radiation pattern of the sonde- 
shell system closely resembles that of an elementary 
oscillating dipole, whose axis makes an angle ¢ with 
the mechanical axis of symmetry (the longitudinal axis) 
of the system. 

Received June 16, 1947. 

* The initial experimental work was done under Section T, 
Division A of the N.D.R.C. 

+ Physicist. 

{t Classified report of the Aberdeen Proving Ground. 

** C. Cranz, Lehrbuch der Ballistiz, dritter band: Experimen- 
telle Baliistik (1927) S. 315-321. 


The modulation envelope for any short interval of 
time is thus: 

(a) Approximately sinusoidal and of a periodicity 
equal to the spin frequency of the shell, w/2z. 

(b) Approximately sinusoidal and of a periodicity 
equal to twice the spin frequency of the shell; that is, 
to 

(c) A’combination of cases (a) and (b) with appro- 
priate relative amplitudes. 

The occurrence of the respective cases (a), (b), and 
(c) depends upon the angle ¢, the polarization of the 
receiving antenna, the trajectory of the shell relative 
to the receiver, and the position of the shell along the 
trajectory. With generalized values of these param- 
eters, case (c) occurs with the relative amplitudes of 
single and double periodicity progressively changing 
during flight. (See section on Receiver and Recorder 
below.) 

Modulation due to varying angle of yaw is super- 
imposed on the spin modulation. This causes a small, 
and relatively low-frequency, periodic variation of the 
apparent instantaneoiis spin frequency. This effect, 
although of course useful for its own sake, is usually 
small and can be eliminated from the determination of 
spin frequency. 

This technique for the measurement of spin of shells 
in flight has had successful use since September, 1941. 

Families of accurate spin curves throughout flight 
for nearly all military shells have been experimentally 
obtained for a wide range of gun elevations and muzzle 
velocities. 

Spin frequency at ¢ = 0 has been used as a convenient 
measure of muzzle velocity. 


THE SPIN SONDE 


During the course of this work various combinations 
of single-ended and push-pull oscillators and two wire 
and loop antennae have been used. Radiated power of 
the order of '/19 watt is sufficient. Higher powers were 
achieved in later developments. 

Typical circuits are shown schematically in Fig. 1a, 
lb, and Ie. 

Vacuum tubes, batteries, and other components are of 
rugged types, developed for the radio proximity fuse. 
The nose of the sond¢ is of insulating material, as lucite. 
Some tests were made with projecting wires (‘‘cat 
whiskers’’) to increase the efficiency of the antenna, but 
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B+ Bt 


LUCITE NOSE 


CIRCUIT & 
BATTERY 


SHELL 


Fic. 2. 


it was found possible to use antennae of both ‘‘two wire”’ 
and loop type lying wholly within the contour of the 
nose. Care was taken to match the weight, moments 
of inertia, and external contour of the standard shell. 
Operation of the oscillator is initiated by ‘“‘setback’’ 
switches that close in the accelerative period in the 
gun barrel. 
ip General arrangement of the sonde-shell combination 
is shown in Fig. 2. 


AXIS OF 
SHELL 


AXIS OF 
EQUIVALENT 
DIPOLE 


Fic. 3. 


“ROTATING SHELL 
TRAJECTORY 


Fic. 4. 


A typical radiation pattern is shown in Fig. 3. By 
careful adjustment ¢ can be made nearly 7/2. Inas- 
much as this is not necessary, it has usually not been 
done. 


THE RECEIVER AND RECORDER 


A Hallicrafter ultra-high-frequency receiver with 
suitable antenna has usually been used for receiving 
and detecting the sonde signal. 

Rotating disc recorders, cathode-ray oscillographs, 
and string oscillographs have been variously used, with 
tuning fork and standard frequency calibrators, to 
make a permanent record of the flight data. 

The amplitude of the radio-frequency voltage in- 
duced in the receiving antenna can be calculated as a 
function of time by the following procedure. (Refer 
to Fig. 4.) 

By neglecting yaw (i.e., assuming the mechanical 
axis of rotation of the shell to lie along the trajectory) 
and using the following notations: 
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LOSS OF SPIN OF PROJECTILES—PART I 


the angle between the electrical axis of the 
equivalent dipole of the sonde-shell sys- 
tem and the mechanical axis 

the angular frequency of rotation of the 
shell 

the distance from the shell to the receiving 
antenna 

right-handed, orthogonal, rotating axes 
with origin at the position of the receiv- 
ing antenna, with the y axis through the 
moving shell, and the x axis lying in a 
plane through the shell axis 

the angle between the x axis and the tra- 


jectory 


the following -amplitude components of the electric 
vector at the receiving antenna are found: 


= 


R 
I 


E, = A/y(cos¢cos a — sin g sin @ cos wt) 
E, = 0 
E, = A/y (sin ¢ sin wt) 


a being a relatively slowly changing function of time 
and A being a constant characterizing the power 
radiated from the sonde-shell system. Ground reflec- 
tion has been neglected here. 

The amplitude of the radio-frequency voltage in- 
duced in the receiving antenna is proportional to the 
numerical value of the projection of this E,, E,, E, vec- 
tor onto its length. 

The most important special case is the one for which 
the receiving antenna lies along the z-axis; then the 
received signal is modulated with time periodicity 
w/a {case (b) of Basis of the Method above], and this 
situation obtains throughout flight. 

Another special case of interest is one with the an- 
tenna along the x axis, a < 7/4, and g = @. Case (a) 
results. 

Fig. 5 exhibits experimental records of cases (a) 
(b), and (c). Determination of shell spin and rate of 


loss of spin can be readily carried out in any one of the 


three cases. 
In using electronic counting, case (b) is the proper 
choice, as indicated above. 


SAMPLE DATA 


Figs. 6a and 6b show some of the early spin-in-flight 
data for a special 57-mm. projectile equipped with two 
different types of antenna. The first was the two-wire 
type, lying wholly within the nose; the second em- 
ployed two hinged dural blades as extensions of the 
antenna. The peak of the trajectory corresponds to the 
point of inflection on these curves. 


APPLICATIONS 


As mentioned above, the radio spin sonde provides a 
method of measurement of muzzle spin and, hence, of 
muzzle velocity. 
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Basic data for the design of spin actuated devices 
in the shell are provided. 

As discussed in Part II by H. P. Hitchcock, skin 
friction coefficients are readily calculated from the 
spin-in-flight data.* 

Fig. 6 illustrates application of the method to the 
determination of the aerodynamic forces on blades ex- 
tending from the nose of the projectile. 

* Suggested by the writer in a private communication to R. H. 
Kent and independently by A. C. Charters and R. H. Kent in an 
Aberdeen Proving Ground report. 
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Part II—Skin Friction Drag 


H. P. HITCHCOCK* 


Ballistic Research Laboratories, Aberdeen Proving Ground 


ABSTRACT 


Analysis of the measured spin of ten types of projectiles leads 
to the evaluation of a skin friction drag coefficient which appears 
to depend on the shape of the projectile as well as Reynolds 
Number. The average value of the skin friction drag coefficient 
is 0.00168. 


DISCUSSION 


I Part I, J. A. Van Allen has explained how the - 
spin of projectiles in flight has been measured by 
means of a radio spin sonde in a dummy fuse. The loss 
of spin Of an experimental 3.3-in. shell had previously 
been measured over a short range by letting a pin pro- 
jecting from the shell body cut notches in a series of 
cardboard screens. The experimental data obtained 
with nine types of high explosive shell equipped with a 
sonde and the one experimental shell with a cutter were 
analyzed by various members of the Ballistic Research 
Laboratories. This analysis led to empirical skin fric- - 
tion drag coefficients, which may be useful in other ap- 
plications besides estimating the loss of spin of projec- 
tiles. 

If w denotes the spin, A the axial moment of inertia, 
t the time, and M, the spin-destroying couple, the spin 
satisfies the differential equation 


A dw ='—Mdt (1) 


R. H. Fowler, E. G. Gallop, C. N. H. Lock, and H. W. 
Richmond have suggested that the spin-destroying 
couple could be represented by the formula 


M, = Kad‘puw (2) 


where K, is called the axial couple coefficient, d is the 
caliber (the diameter of the bore and approximately 
that of the projectile), p the air density, and u the veloc- 
ity of the projectile relative to the air. R. H. Kent ap- 


* Ballistician, Ext. Ball. Lab. 


plied this formula to obtain the solution of Eq. (1): 
= exp[(—d*/A) fo'Kapu dt] (3) 


if wo is the spin when ¢ = 0. 

The spin may be measured at two or more times 
along a trajectory, as explained before. The caliber of 
the gun is known. The axial moment of inertia of the 
shell can be determined by swinging it about its axis on 
a torsion pendulum and comparing its period with that 
of two different cylinders swung in the same way. 
The air density and velocity along a trajectory can be 
computed by numerical integration for a given muzzle 
velocity, angle of departure, and ballistic coefficient.' 
Hence, the axial couple coefficient can be computed by 
Eq. (3), providing this coefficient is considered con- 
stant. 

A. C. Charters and R. H. Kent showed that the skin 
friction drag may be approximated by the relation 


Dy = 4uM,/d*w (4) 
and that it may also be expressed as 
Dy = Cor’ (5) 


where Cpr’ is a skin friction drag coefficient and S’ the 


surface of the projectile exclusive of the base. It fol- 
lows from Eqs. (2), (4), and (5) that 
Cor’ = 4K4d?/S’ (6) 


Since the surface of the projectile can be calculated from 
its dimensions, the skin friction drag coefficient can be 
calculated from the axial couple coefficient. 

According to von Karman and Prandtl,’ the skin 
friction drag coefficient pertaining to the motion of a 
turbulent boundary layer on a flat plate approximately 
satisfies the empirical relation 


Cor’ = 0.037R-°-? (7) 


38 
i 
3.3 
—— 
90 
i 
ay 
1 
ae 
4 
mu 
in 
air 
pet 
an 
4 un 
die 
th 
pr 


pursuing 
appreci- 


e 4.7-in, 
H. Sreb, 


prepara. 
the firing 


1): 
(3) 


times 
iber of 
of the 
xis on 
h that 

way. 
an be 
1uzzle 
ient.! 
by 
con- 


skin 


(4) 


(5) 


the 
fol- 


LOSS OF SPIN OF PROJECTILES—PART II 


TABLE 1 
Axial Couple and Skin Friction Drag Coefficients I 
Skin 
: Reynolds Axial Friction 
Muzzle Time Number Couple Drag 
Velocity, Elevation, Interval, 0.45 WLp/e, Coefficient, Coefficient, 
Shell Ft. per Sec. Deg. Sec. logioR Ka Cor’ 
75-mm. HE M4s 2,014 10.8 1.1-18 6.6057 0.00587 0.00175 
3-in. HE M42 2,747 60 6.2-36.5 6.2916 
2,718 45 5.7-41.6 6.3540 
2,711 30 1.5-29.1 6.4394 
2,722 15 0.3-21.0 6.6970 
Avg. 6.4455 0.00585 0.00218 
3.3-in. Exp. T68 2,015 —0.4 0.02—-0.51 6.8706 0.0048 0.00162 
90-mm. HE M71 2,736 60 9.2-53.2 6.2957 
2,739 45 5.8-45.2 6.4124 
2,737 30 1.7-9.7 6.7655 
2,737 15 2.7-24.3 6.6669 
. Avg. 6. 5351 0.0059 0.00198 
4.5-in. HE M65 2,379 15 0.8-28.4 6.8786 
2,380 45 1.2-68.9 6.6154 
1,870 20 1.4-27.4 6.7897 
Avg. 6.7612 0.00535 0.00163 
4.7-in. HE M73 3,010 9.25 5-23 6.9344 0.00507 0.00165 
approx. 30 5-20 6.8490 0.00447 0.00147 
30 20-40 6.6197 0.00587 0.00193 
30 40-60 6.7040 0.00528 0.00173 
80 5-20 6.6869 0.00405 0.00133 
80 20-40 6.1741 0.00602 0.00198 
80 40-60 5.7362 0.00696 0.00229 
80 60-80 5.9575 0.00487 0.00160 
80 80-100 6.4243 0.00401 0.00132 
80 100-120 6.7248 0.00348 0.00114 
Avg. 6.4811 0.0050 0.00164 
5-in. AA Mk31 2,600 10 9.0-20.1 6.7394 0.00494 0.00180 
approx. 20 15.8-35.0 6.6418 0.00567 0.00207 
45 18.0-40.0 6.4099 0.00567 0.00207 
60 21.0-79.5 6.3768 0.00477 0.00174 
* Weighted Avg. 6.4955 0.00515 0.00188 
155-mm. HE M101 2,114 18 6.3-27.5 6.8524 
2,813 13 6.9-26.4 6. 8602 
Avg. 6.8563 0.00395 0.00147 
8-in. HE M106 1,948 18 0-24 7.0277 0.0045 0.00172 
240-mm.HE M114 2,310 17 7.0-10.2 7.1926 
1,487 20 3.3-27.8 7.0251 
1,452 20 3.5-26.7 7.0251 
2,292 20 7.7-34.4 7.0631 
2,316 20 6.8-31.6 7.0634 
Avg. 7.0739 0.0026 0.0097 


where R is Reynolds Number defined by the for- 
mula 


R = WLp/c (8) 
in which W is the wind, L the length of the plate, p the 


air density, and o the air viscosity (a function of tem- 
perature). In dealing with projectiles, if we take 


W = (uw? + (9) 


and L as the length of the projectile, a factor less than 
unity should be inserted to account for the variation in 
diameter, since a shell is not cylindrical. The factor 
that makes the formula agree with the average of the 
present results is 0.45. 


The computed axial couple coefficients and skin fric- 
tion drag coefficients are listed in the table, which also 
gives the projectiles fired in the tests, the muzzle 
velocity, the elevation of the gun, the shortest and 
longest time of flight for which the spin was measured, 
and the average value of Reynolds Number over the 
part of the trajectory considered. From two to ten 
shots were fired from each gun. The long trajectories 
of the 4.7-in. HE shell at 30° and 80° were divided into 
sections as indicated. If the results for each of the ten 
projectiles are assumed to have equal weight, the mean 
value of the common logarithm of Reynolds Number 


‘ is 6.7153 and the mean value of the skin friction drag 


coefficient is 0.00168. 
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On the graph (Fig. 7), the computed values may be 
compared to the curve determined by von Karman’s 
formula (7) for a flat plate. Although the skin friction 
drag coefficient seems to depend on Reynolds Number, 
this relation appears to be affected by the shape of the 
projectile. 
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Letters to 


Dear Sir: 

In his article ‘The Ring Airfoil in Nonaxial Flow” in the Sep- 
tember, 1947, issue of the JouRNAL OF THE AERONAUTICAL ScI- 
ENCES (Vol. 14, No. 9, pp. 529-530), Mr. Ribner found that the 
lift of a ring airfoil was given by [Eq. (5)]: 

L = 1/[(1 + (e/4R)] p V2x*Rea 

It is interesting to note that practically the same result had 
been obtained by the writer by considering the simple limiting 
case of a ring airfoil as R > ~, while c remains small. 

In this case (following the notation given in Fig. 1 of Mr. Rib- 
ner’s article), the local angle of attack for each section of the upper 
half of the ring may be written as (a sin @), and the mutual inter- 
ference between adjacent sections becomes negligible as R > ~, 
Then 


1 dC 
L—~2f" sin? 6 


where (dC;/da) may be taken as the constant lift curve slope of 
the profile. Consequently, 
L = V*cR(dC,/da)ar - ‘ 
If one takes (dC,/da) = 2x, this expression agrees with Rib- 
ner’s Eq. (5) for the limiting case where (c/R) — 0. 
If one takes 


(dC, /da) = 2x[A’/(A’ + 2)] 
this expression also agrees with Ribner’s Eq. (5) for the same as- 
pect ratio he defines as A’ = 8R/zc. 
E. V. LAITONE 
Associate Professor 
University of California 


the Editor 


Dear Sir: 

In the opening paragraph of his paper on the “Effect of Wind 
Gradient on Glide and Climb”’ in the June issue of the JouRNAL, 
Bernard Etkin remarks that “incorrect conclusions’’ on the sub- 
ject were contained in an article written by the undersigned for 
Canadian Aviation in October, 1943. 

This writer can lay no claim to originality with respect to 
these conclusions, since they have been published on at least 
two other earlier occasions. In the 1936 edition of Air Publica- 
tion No. 129, the official flying manual of the Royal Air Force, 
the case is stated and illustrated clearly in the chapter on flying 
technique. 

Dr. M. A. Garbell included the question in his Basic Aeronau- 
tical Meteorology written in 1941 as a textbook for the Boeing 
School of Aeronaaitics. Dr. Garbell writes of flight through a 
wind gradient. ‘An airplane makes dynamic gains while climb- 
ing against the wind and experiences dynamic losses while climb- 
ing with the wind. Similar gains are made while gliding with the 
wind and losses while gliding against the wind.” 

Analyzing the problem from energy considerations aione, it is 
evident that, if the wind velocity decreases more rapidly than the 
gliding airplane can accelerate along the glide path, then, to 
maintain his air speed, the pilot must obtain energy from the 
only possible source in a glide—namely, by losing additional al- 
titude. 

Finally, experienced pilots can testify that, while approaching 
a landing field through a steep wind gradient, the airplane sinks 
rapidly during the last 200 ft. of descent. 


Luioyp H. SLOAN 
Piasecki Helicopter Corporation 
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ABSTRACT 
Uni- 

This paper describes the derivation of an optimum inlet Mach 

III, Diy Number for the first stages of an axial-flow compressor based on 
_’ | free vortex design. The optimum Mach Number is related to 
T, Berlin, the given inlet conditions, dimensions, and tip speed and specifies 
the maximum amount of work that can be done on the fluid and 
the corresponding free vortex velocity diagrams at any radius 

in a stage where Mach Number effects are critical. 

A graphical method is described for simplifying the construc- 
tion of free vortex velocity triangles for a stage at any radius of 
flow. This method is derived from the definition of free vortex 
flow. 

It is shown that it is possible to design several successive stages 
of Wind using the same rotor blade form and the same stator blade form 
OURNAL, by merely changing the spacing and root stagger angle. This 
the sub- may be done without compromising the desired load distribution 
sned for — among stages and without violating the laws of free vortex flow. 
List OF SYMBOLS 
at least 
ublica- a = acoustic velocity, ft. per sec. 

Force, a = acoustic velocity at stagnation conditions, ft. per sec. 
1 flying Cm = axial component of velocity, ft. per sec. 

C, = profile pressure coefficient dimensionless = 

ronau- (p — 
Boeing (Cp)o = profile pressure coefficient at zero Mach Number— 
ugh a dimensionless 
climb- (Cp)m = profile pressure coefficient at Mach Number = M— 
climb- dimensionless 
th the Cul = tangential component of ¢, ft. per sec. 
¢y = tangential component of co, ft. per sec. 
, it is Acy = vector change in velocity, ft. per sec. 
in the ra = absolute velocity entering stage (entering rotor), ft. 
n, to per sec. 
n the Ce = absolute velocity leaving rotor, ft. per sec. 
al al- C3 = absolute velocity leaving stage (leaving stator), ft. 
per sec. 
ching D diameter of flow section, ft. 
sinks D) = diameter of flow section where velocity diagram is 
symmetric, ft. 
D, = inside diameter of casing, ft. 
- g = acceleration of gravity, 32.2 ft. per sec.? 
tion H = isentropic total enthalpy rise of one stage, ft.lbs. per 
Ib. 

K_ = radial parameter—dimensionless = D/Do 

M, = Mach Number of flow entering stage—dimensionless 

M. = critical Mach Number—dimensionless 

= revolutions per minute 

p = static pressure, lbs. per sq.ft. 

pi: = static pressure in front of cascade, lbs. per sq.ft. 

Q = flow volume entering stage, ft.* per sec. 

r = radius of flow section, ft. 

U- = peripheral velocity of rotor blade, ft. per sec. 

Presented at the National Aircraft Propulsion Meeting, I.A.S., 
Cleveland, Ohio, March 28, 1947. 

* Division Engineer, Aerodynamics Division, Research and 
Development Department. 
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A New Method for the Aerodynamic Design 
of Multistage Axial-Flow Compressors 


HUNT DAVIS* 
Elliott Company 


U, = peripheral velocity of rotor blade at diameter D, ft. 
per sec. 

U,. = peripheral velocity of rotor blade at outside diameter, 
ft. per sec. 

W = work done on fluid by a rotor. row, per Ib. of air, 


ft.lbs. per Ib. 
relative velocity at entrance to rotor row, ft. per sec. 


= 

We = relative velocity leaving rotor row, ft. per sec. 

a = absolute angle of fluid entering rotor row, degrees. 
See Fig. 1 

a, = absolute angle of fluid leaving rotor row, degrees 
See Fig. 1 

a; = absolute angle of fluid leaving stator row, degrees. 
See Fig. 1 

Aa = deflection angle of fluid in stator, degrees. See Fig. 1 

8: = relative angle of fluid entering rotor row, degrees. 
See Fig. 1 

B: = relative angle of fluid leaving rotor row, degrees. See 
Fig. 1 

Ag = relative deflection angle of fluid in rotor row, degrees. 
See Fig. 1 

n = isentropic stage efficiency—dimensionless 

v = ratio of root diameter to tip diameter f— limensionless 

p = density of fluid, Ibs. sec.* per ft.‘ 

= flow coefficient = ¢»/Us—dimensionless 

v = pressure coefficient = 2gH/U»*—dimensionless 

¥/n = work coefficient—dimensionless 


INTRODUCTION 


few RAPID RISE of the gas turbine to prominence 
during the war years was accelerated because of its 
suitability for power in jet-propelled aircraft and guided 
missiles. However, the real advent of the gas turbine 
might have come years earlier but for two factors. One 
of these was the lack of suitable high-temperature al- 
loys, which has been only recently overcome by the re- 
markable progress in high-temperature metallurgy. 
The other was the need for improvement in efficiency 
and reliability of high-performance compressors. 

There is no doubt that we are still very much in the 
dark about the aerodynamic intricacies of centrifugal 
and axial-flow compressors. This paper will be de- 
voted to a few aspects of axial-flow compressor design. 

For lack of space, we are forced to neglect entirely a 
presentation of the effects of Reynolds Number, aspect 
ratio, leakage, and secondary flow effects. We will 
deal here with an inlet Mach Number discussion and a 
graphical method for making calculations of blade set- 
tings. It will be shown that it is possible to use the 


+ In this paper “tip” is equivalent to “outside diameter of flow 
path” and “root” is equivalent to “inside diameter of flow 
path.” This refers to both rotor and stator rows. 
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Fic. 1. Typical velocity vector diagrams for one axial-flow 
compressor stage. 


same blade design for successive stages even though the 
stage loadings are not constant. 


VELocITY VECTOR DIAGRAMS 


Let us consider the velocity vector diagrams for an 
axial-flow compressor stage. It is assumed that all the 
streamlines of the flow lie on cylinders concentric with 
the axis of rotation. On such a cylindrical surface the 
velocity vector diagrams shown in Fig. 1 may be ob- 
served. 

It can be shown that for frictionless irrotational flow 
through a stage, the axial component, c,, does not 
vary. along the radius. (By axial component, c,, is 
meant the average value that ignores periodic variations 
due to blade wakes.) In Fig. 1, c, has been assumed 
constant through one stage as well. This is very nearly 
the actual case for one stage since the usual decrease 
in annulus area nearly compensates the increase in 
fluid density. We define 


Ac, = Cu2 nara Cul (1) 


Recalling that the work W, per pound of fluid, done by 
the rotor element is 


W = Ude,/g (2) 
and that the isentropic stage efficiency is 
n = H/W (3) 


we may prove that 
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¥/n = 2U Ac,/ Ue? (4) 
This follows from the definition of y which is 
y = 2gH/U? (5) 


the pressure coefficient. Up» is the peripheral blade 
speed at a reference diameter. y/n is a dimensionless 
group, called the “work coefficient,” relating the change 
in tangential velocity to the velocity of the blade. 

We may define a ‘‘flow coefficient,’”’ y, as the ratio of 
the axial component to the velocity of the blade at a 
reference diameter. 


Cm/ Uo (6) 


¢ will be constant for the whole radial height of the 
blade, since ¢,, is constant. Similarly, ¥/y is constant 
over the blade height, because the work done per 
pound of fluid by a blade element is to be constant all 
along the blade. 


FREE VORTEX THEORY 


In order to design a highly efficient compressor, con- 
sideration must be given to the various types of blading 


- that may be used. The flow through a compressor is, 


spiral in general, since at every point the absolute veloc- 
ity consists of an axial component and a tangential 
component. (The radial component may be caused to 
be negligibly small.) The radial variation of these 
components determines the type of blading as well as 
the radial distribution of pressure. The radial distribu- 
tion of pressure is related to the tangential velocity 
component c, for perfect incompressible fluids by 


in which p is the pressure, p is the density, and 7 is the 
radius. The last term on the right-hand side of Eq. 
(7) may be dropped, since it is always small compared 
to c,”/r for axial flow compressors with free vortex de- 
sign. 

If the flow entering a blade row is of the potential 
type and if the blade is shaped so that the tangential 
component at the exit varies with the radius according 
to some law, there may be potential flow at the dis- 
charge of the blade row, provided that the blades are 
designed to deflect the fluid in both the radial and tan- 
gential directions. The amount of the radial deflection 
depends on the prescribed radial variation of the tan- 
gential component and upon the variation in density 
due to the compressibility of the fluid. Only when 
rc, = constant is no radial deflection required for in- 
compressible fluids and is the axial velocity constant. 
(When rc, = constant and the fluid is compressible, 
then there will be a radial deflection to satisfy the con- 
tinuity condition with the radial variation in density. 
However, c,, will remain constant.) This law forms 
the definition of ‘‘free vortex’” flow for axial-flow com- 
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pressors with either compressible or incompressible 
fluids. 


1Cy 


constant in each radial wre (8) 
Cm 


constant 


Of the many types of blading that could have been 
chosen, the free vortex type has been chosen for this 
analysis for its simplicity and because it satisfies the 
potential flow problem for a perfect incompressible 
fluid without radial cross flows. The fotegoing remarks 
are not an attempt to justify free vortex design but 
rather serve as an introductory review for the methods 
described below. 


Macy NuMBER LIMITATIONS AND DESIGNING FOR THE 
OptimuM INLET MacH NUMBER WITH FREE VORTEX 
FLow 


The characteristic thinking of aircraft-engine de- 
signers along the lines of ‘“‘more power in less space” 
when applied to a compressor leads to the conclusion 
that the rotor tip speed should be as high as is feasible. 
Two things limit the maximum tip speed of a compres- 
sor: the high stresses resulting from centrifugal force, 
and prohibitive shock losses that arise when the veloc- 
ity of the fluid relative to a blade approaches the acous- 
tic velocity. Accordingly, each stage must be designed 
with the above effects in mind. 

In order to minimize shock losses in a stage, an opti- 
mum concept has been developed which will aid 
in designing the stages where the Mach Number is 
high. 

Whenever the local velocity along any streamline 
approaches the acoustic velocity at the local tempera- 
ture, the flow deviates from the intended pattern and 
a localized loss in efficiency is realized. 

The Mach Number at any point on a streamline is 
M = c/a, where c is the flow velocity at the point and 
a is the acoustic velocity at the point. 

The discussion of Mach Number effects will be di- 
vided under two headings. 

(1) Calculation of the critical Mach Number for 
any profile from the entrance conditions. 

(2) Calculation of an optimum design for the first 
stages of a machine where the Mach Number limita- 
tions completely determine the performance of a 
stage. 

In order to calculate the local Mach Number of a flow 
through a row of blades, a relationship is necessary be- 
tween the maximum local velocity and the undisturbed 
entrance velocity for a given cascade with a given en- 
trance angle. 

It is obvious that the maximum local Mach Number 
is greater than the Mach Number of the undisturbed 
stream. The maximum local Mach Number occurs at 
the place on the profile where the local pressure is a 
minimum. 
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Fic. 2. Profile pressure distribution around an airfoil in cascade 
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Fic. 3. Relation between the profile pressure coefficient and 
critical Mach Number.! 


A typical pressure distribution curve for an airfoil 
in cascade is shown in Fig. 2, where the pressure coef- 
ficient C, is plotted against the position along the 
chord. C, is defined by 


C, = (b — (9) 


where /; is the pressure in the undisturbed flow up- 
stream of the cascade. The point of maximum velocity 
is found where the value of C, is a minimum. For a 
given single airfoil, with a given angle of attack, C, 
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Fic. 5. Relation between the velocity ratio and Mach Number 
for air.? 


8 3 


changes with the Mach Number of the upstream flow 
according to an equation given by von Karman! 


(C,)o 


\v (; + 
(10) 


(Cy) m 


in which (C,)m is the coefficient to be found, and (C,)o 
is the coefficient at Mach Number zero. ; 

The critical Mach Number is defined to be the Mach 
Number of the upstream fluid when the maximum local 
Using Eq. 10 and one relating 


Mach Number is one. 


pressure ratio with Mach Number changes, von Karman 
obtains the curve of Fig. 3.1 This curve relates the 
critical Mach Number for a single airfoil to the mini- 
mum C, value determined from tests at low Mach Num- 
bers. Although the derivations are for single airfoils, 
they are used here for cascades. The discrepancies are 
of small but uncertain magnitude and are probably 
mainly due to the overall change in main stream veloc- 
ity which is produced by cascades and not by single 
airfoils. In any, case, certain assumptions have been 
‘made in the derivations, and different writers have 
found other functions than the one shown by Fig. 3. 
von K4rman’s curve was taken because it seems to be 
the most conservative of all the ones that could be 
chosen and fits the single airfoil data better. 

It is recommended that machines shall be designed 
so that the inlet Mach Number for a stage is never 
greater than the critical Mach Number according to the 
curve of Fig. 3. 

In the free vortex design of the first stages of a ma- 
chine with a given flow, the maximum work that can 
be put into compressing the air usually depends en- 
tirely on the stage dimensions and the allowable Mach 
Number of the blading. 

Referring to Fig. 4, which shows characteristic veloc- 
ity diagrams at three different radii of one stage, we 
see that for the rotor blades the maximum inlet Mach 
Number occurs at the tip, while for the stator blades it 
occurs at the root. We shall assume that an optimum 
first stage design is one where 


(11) 


This will make the two local maximum Mach Numbers 
equal if: (1) the temperatures at the stagnation points 
relative to the respective blades are equal at the rotor 
tip and stator root; (2) the values of (Cy) min. are the 
same at rotor tip and stator root. For the time being 
we shall assume that these conditions are fulfilled. 

The expression for (w;) ;, has been derived and is 


4 1 2 
(W)tip = (€)root = {3 


gw 2 
ze 
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The corresponding Mach Number M is found by 
dividing w, by a and then reading M from the curve of 
Fig. 5.? 

It can be shown with Eq. (12) that for every Q, 
D,, U, and W value there is an optimum finite v for 
which w, is a minimum. This optimum yr value varies 
from about 0.65 to 0.80 for various compressors. The 
optimum » is obtained by differentiation of w,; with re- 
spect to v. 

It should be pointed out that the usefulness of the 
optimum concept does not lie in minimizing w, for a 
given amount of work, but-in obtaining the maximum 
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stage work for a given allowable maximum value of 
and for a given flow. 

Using the condition with optimum v concept we can 
construct diagrams, similar to Figs. 6a ani 6b, for each 
given maximum allowable w, value. Figs. 6a and 6b 
have (w:),ip = 650 ft. per sec. Fig. 6a shows the op- 
timum diameter ratio v as a function of Q/D,. Fig. 
6b shows the relationship between the maximum work 
W, the tip speed U;,, the optimum »v, and the parameter 

Fig. 6b involves the quantities Q, N, W, v and D,. 
When N and any two other quantities are specified, 
the remaining ones are determined. This defines the 
complete first stage design, which is optimum for an 
inlet maximum w, = 650 ft. per sec. 

Upon reviewing the two preliminary assumptions, 
we see that the first is not true because the air entering 
the stator row has a different stagnation temperature 
from the air entering the rotor row. This difference is 
due to the work done by the rotor and to the difference 
in reference frames. 

The validity of the second assumption depends on 
the profiles that are used at the root and tip sections 
of the stator and rotor, respectively. 

The fact that these assumptions are not met requires 
that the inlet Mach Number be slightly less than the 


volved. 

When a stage-by-stage design calculation is carried 
out, the Mach Number must be watched carefully at 
the tip of the rotor blades and the root of the stator 
blades for several stages‘after the first until it has passed 
from the critical range. 

If the Mach Number at the rotor tip is too high, it 
may be lowered by moving the symmetrical section 
along the blade toward the tip. This increases the 
Mach Number at the stator root, however, and it must 
be rechecked. A lower value of the maximum Mach 
Number should be used if there will be air temperature 
changes at the machine inlet or other uncertainties. 
In general, the position of the symmetrical section in the 
first stages is determined entirely by Mach Number 
considerations. After the first stages, it should be as 
close to the root of the blades as possible without the 
highest inlet Mach Number exceeding the critical 
value at any point on the blade length. This is de- 
sirable because the stalling characteristics of rotor and 
stator blades determine approximately a maximum al- 
lowable turning angle for a given profile cascade. Fora 
given turning angle of a rotor section, the amount of 
work that can be done in the stage will increase if the 
symmetrical section is moved toward the center of the 
machine. However, if the K value at the root is made 
greater than 1.0, the stator tip deflections usually be- 
come excessive. Usually K = 1.0 at the root is opti- 
mum for maximizing the stage work with turning angle 


limitations. 
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critical value to allow for the small uncertainties in- 


80 
w.=650 FIL/SEC. 


= \ 
= 68 
a 
64 
77 6 8 1.0 12 


CFS/IN? 


Fic. 6a. Relationship between the optimum » value and Q/D,? 
for free vortex flow with w, = 650 ft. per sec. 
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Fic. 6b. Curves giving design values when inlet Mach Number 
is optimum. [(wi)tip = (C2)rooe = 650 ft. per sec. 


VELOcITY DIAGRAMS FOR FREE VORTEX BLADING 


Velocity vector diagrams for free vortex flow for a 
stage are shown in Fig. 4. The upper diagram shows 
the velocities for a radial element where the vectors are 
symmetrically placed. This is known as the “mean” 
or ‘‘symmetric”’ section, and it is the section where the 
degree of reaction is 50 per cent. The diameter of the 
mean section is Dy. Other diameters are related to Do 
by the equation D = KD). The lower two diagrams of 
Fig. 4 represent the vectors for diameters where K = 
1.2 and 0.8, which are characteristics of tip and root 
conditions, respectively. These three diagrams fulfill 
the conditions set forth for free vortex flow, which are: 
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Fic. 7. Diagram showing use of velocity vector work sheet. 
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Fic. 8. Rotor deflections as a function of ¢/y and radius for 
¢g = 0.40 


(1) Uc, is constant along the radius, (2) U Ac, is con- 
stant along the radius, and (3) c,, is constant along the 
radius. 


The assumption is made that c,, is constant through 
the stage, as well. 

Certain quantitative relationships are obvious from 
Fig. 4." For instance, the variation of rotor deflection 
A@ is larger from root to tip than the variation of stator 
deflection Aa. Also, the variation in relative velocities 
W}, W2 is greater than the variation in absolute velocities 
Ci, C2. 

A useful diagram for aid in calculating vector quan- 
tities is the ‘‘Velocity Vector Working Sheet’? shown 
in Fig. 7. Its construction follows directly from the con- 
ditions of free vortex flow. The w/y lines are lines of 
constant “‘loading,” and every c/U,) and vector 
must end on the same y/y line to satisfy the condition 
U Ac, = constant. 


CALCULATION METHOD FOR REPEATED USE OF THE 
SAME BLADING IN SUCCESSIVE STAGES 


Frequently, a compressor is designed in such a way 
that the blade form in every stage is different from all 
the other stages, thus requiring a different die for each 
blade row. It would be advantageous economically to 
be able to use the same blade design in several stages 
without compromising the desired load distribution be- 
tween stages and without violating the laws of free 
vortex flow. It can be shown that such a design is 
possible. We will proceed to describe the methods used 
to do so, assuming that the only difference in blading 
between several successive stages is a reduction in length 
obtained by cutting off the free end and a different angle 
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setting of the blade in the root fastening. This will 
be done with both rotor and stator blades. 

We will need charts such as those shown in Figs. 8 
and 9. These figures show the variation in turning 
angles AB and Aa with K, the radius parameter, and 
y/n, the work coefficient for one value of ¢, the flow 
coefficient. These curves and ones for other ¢ values 
can be calculated from the geometry of the vector 
triangles and the derivations are omitted here. 

Let us assume that by some means, such as the opti- 
mum inlet Mach Number theory, we have arrived at 
the rotor and stator blade design for the first stage of a 
multistage machine. We shall also assume that the 
diameters of casing and hub have been chosen for each 
stage. We wish to design the second stage, using the 
same blades as the first with different fastening set- 
tings. 

A’ paper template is made and is used as shown in 
Fig. 10 (assuming, of course, that the ¢ value of the figure 
shown is correct). The template is cut to the contour of 
the y//n line for the first-stage performance. A vertical 
index line is marked corresponding to the K,,o¢ value. 
This template is identified with the rotor blade of the 
first stage. If the same rotor blade is to be used for the 
second stage, then this template may be used with the 
appropriate curve sheet for the-new ¢ of the second 
stage. Moving the template horizontally across the 
curve sheet corresponds to moving the symmetrical 
section in and out along the blade length. Moving the 
template up and down by small amounts on the curve 
sheet corresponds to changing the absolute stagger 
setting of the blade by small amounts. The second- 
stage work coefficient, ¥/n, should be higher than that 
in the first stage to the extent permitted by Mach Num- 
ber considerations. Hence, the template is adjusted 
so that its contour line fits some satisfactory y/n line, 
the index line, of course, remaining vertical. The index 
line then indicates a new K,,., value for the second-stage 
design. From this quantity and the indicated ¢ value, 
the second-stage rotor design may be laid out in full on 
the velocity vector sheet. Although the exact ¢ value 
depends on A,oot, the template manipulation process is 
not at all sensitive to changes in gy. Hence, a ¢ value 
may be assumed to start in picking the appropriate 
curve sheet on which to use the template. 

The second-stage rotor -design prescribes a certain 
Ca value at every radius. The first-stage stator design 
prescribed a certain c,, at its discharge. In general, 
the two values will not be matched. A correction must 
be made to match the stator discharge to the rotor inlet 
of the following stage. The adjustment is made by 
making a small,change in the absolute stator stagger 
setting so that the stator blades produce the required 
deflection corresponding to the prescribed rotor inlet 
Ca value. This correction does not upset “free vortex” 
conditions but only changes the velocity heads slightly. 
These computations are simplified by the use of the 
velocity vector working sheets. For the remaining 
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Fic, 10. Diagram PR use of template on Fig. 8 to facilitate 
design using same blading in successive stages. 


stage calculations, the most convenient sequence of 
computation is to calculate both rotor rows on either 
side of a stator row before calculating the stator row. 
This immediately determines the necessary stator de- 
flections and stagger settings. 

Although c,, and the flow area change continuously 
through the machine, we assume they are constant at 
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an average value with small discontinuities between 
stages. This simplifies calculations and works well. 

In matching one stator row to the following rotor 
row, we must consider the change in absolute value of 
the “root,” ‘‘middle,”’ and ‘‘tip’’ dimensions when using 

- the velocity vector working sheets. A useful relation is 


(CurDoK root) stage (CurDoK root) stage B (13) 


After a few stages after the first have been computed, 
the template will no longer fit any ¥/n curve that is as 
high as is needed for higher stage loadings. Hence, a 
new blade and corresponding template are designed. In 
general, the template procedure is not necessary for 
stator blades because of the flatness of the curves of 
Fig. 9. The need for redesign of stator blades in going 
from one group of stages to the next will be self-evident 
from the velocity vector work sheets. Particularly, 
this need is seen when the rotor blade design is changed 
between two stages. 

Care must be taken to see that the rotor deflections 
at the root do not exceed the stalling deflections for 
the camber and solidity values used. At the stator 
tips, the deflection should be definitely below the stall- 
ing value, as indicated by cascade tests. A reasonable 
value is that the deflection at the stator tips should not 
exceed 0.7 of the stalling deflection as measured by a 
stationary cascade test. The stalling situation is re- 
lieved at the rotor roots and aggravated at the stator 


tips by the radial displacement of the profile boundary r 
layer along the blade contours. This displacement of 
low relative velocity fluid by centrifugal force causes 
an accumulation of low-energy fluid at the tips anda 
replacement at the roots of high-energy fluid. The low- 
energy fluid has incorrect flow angles entering the stator 
blade row, and, consequently, stalling is aggravated. 


CONCLUSION 


It is hoped that the reader does not think that we |. - 
have jumped into the middle of the problem of com- ont 
pressor design without due attention being focused upon J jhe sp 
the preliminary design calculations. The author has | shown 


assumed that the obvious problems of the determina equati 


tion of the number of stages, the selection of blade pro- | *4 . 
files, the calculation of stage and overall efficiency, the ae 
variation of root and tip diameters from stage. to stage, 
the variation of stage loading, and other such problems As | 


have already been solved by other means. The author § Savin 


plans to discuss some of these other points in a future hn 
stre 

paper. 
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SUMMARY 


Starting with the Navier-Stokes equations of motion for vis- 
cous fluids, Prandtl’s argument regarding the orders of magni- 
tudes of quantities in a laminary boundary layer is repeated for 
the special case of a yawed cylinder of infinite length. It is 
shown that for this case the chordwise flow is given by the same 
equations as for unyawed (plane) flow about the same cylinder 
and that the spanwise boundary-layer flow can be calculated by 
integration of a linear second-order differential equation. This 
result leads to useful conclusions—e.g., that the location of 
laminar separation of the chordwise flow is independent of yaw. 

As an example, the flow over a flat plate, at zero incidence, 
having a sweptback leading edge, is treated. It is found that 
the boundary-layer flow is always in the direction of the external 
stream and has the Blasius velocity profile—i.e., it is unaffected 
by sweepback of the leading edge. As a second example, a cer- 
tain class of cylinders suggested by Prandtl is considered in the 
yawed condition. The flow about these in the plane case is 
known, the spanwise boundary-layer flow is calculated approxi- 
mately, and some characteristic features are presented graph- 
ically. 


INTRODUCTION 


— EXPERIENCE with sweptback wings has 
revealed some aerodynamic phenomena that are 
not readily explained by theories based on the friction- 
less-fluid approximation. These phenomena are re- 
flected in nonlinear and irregular behavior of force and 
moment curves of sweptback or yawed wings, espe- 
cially at large angles of attack, and are apparently as- 
sociated with profound effects of viscosity in the bound- 
ary layers of the wings. It seems desirable to begin a 
study of the boundary layers of finite-span wings, es- 
pecially yawed and sweptback, and the effects of three- 
dimensional boundary-layer flows on the characteris- 
tics of such wings. It is also of interest to study such 
boundary layers in tle hope of determining the effects 
of yaw and sweepback on the important boundary- 
layer phenomena of transition and separation. 

As a first step in this program, the laminar boundary 
layer of an infinite yawed cylinder is considered in some 
detail here. It will be shown that the simplifying as- 
sumptions made here—i.e., laminar flow over a body 
of uniform section, infinite in extent—permit an essen- 
tial simplification of the three-dimensional boundary- 
layer equations, so that some problems of considerable 
interest can be solved. 

In the case of frictionless fluid the flow about yawed 
infinite cylinders has been constructed! by considering 
separately the effects of the stream components per- 
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pendicular and parallel to the axis of the cylinder. 
These flows can be superimposed; in fact, the axial 
component is uniform throughout the region and is 
equivalent to a translation of coordinate axes. For a 
viscous flow, of course, this extremely simple device 
cannot be used; nevertheless, it will be shown that a 
simplification results which permits two-dimensional 
(unyawed) flow solutions to be employed in the calcula- 
tion. 


I. EQUATIONS 


The boundary-layer equations for a yawed infinite 
cylinder are obtained from the Navier-Stokes differen- 
tial equations of motion by repeating, with the neces- 
sary alterations, L. Prandtl’s argument concerning the 
magnitudes of the various terms for small viscosity 
(large Reynolds Number). Following Goldstein,? one 
introduces curvilinear coordinates x and z, tangent and, 
perpendicular, respectively, to the surface of the cylin- 
der in planes normal to the cylindrical axis. Thus, x 
and z are suitable coordinates to describe the plane 
flow in the case of the unyawed cylinder. A third 
coordinate y is then introduced in the direction of the 
cylindrical axis (see Fig. 1). One defines the corre- 
sponding velocity components in the directions x, y, 
andzasu,v,andw. ° 

An argument, analogous to Prandtl’s,?"* can now be 
constructed, based on the observation that uw and v must 
vanish at z = 0 and must attain their potential-flow 
values u; and v, at the outer edge of the boundary 
layer, whose thickness 6 is small compared with the 
pertinent dimensions of the cylinder. Thus, 0u/dz 
and 0v/dz must be of order 07u/0z? and 0*v/dz? 
must be O(5-?), etc.; while u, v, 0u/Ot, dv/Ot, Ou/odx, 
dv/dx, etc., must be O(1), in the boundary layer.t 


denotes the time. 


Fic. 1. Sketch showing coordinate system used. 
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Moreover, the cylinder is uniform and infinite in ex- 
tent, so that none of the physical quantities that de- 
scribe the flow can be functions of y. Thus the equa- 
tion of continuity becomes: 


(Ou/Ox) + (Ow/dz) = 0 (1) 


and all derivatives with respect to y in the equations of 
motion vanish. Eq. (1) indicates that Ow/dz is O(1); 
hence, w is O(6), as are Ow/Ot, Ow/Ox, 0?w/dx?, etc. 
The Navier-Stokes equations, under these special 
conditions and with terms of order 6 being neglected, 
are reduced to the following, where v denotes the kine- 
matic viscosity and p the density: 


ou Ou Ou 1 Op 
ov Ov Ov 077 
(3) 
dp/dz = 0 (4) 


It is seen that Eqs. (1), (2), and (4) do not involve 
either the velocity component v or the coordinate y; 
they are, in fact, identically the same as for the two- 
dimensional (plane) case. They can therefore be 
solved, for given boundary conditions, for u = u(x, z, t) 
‘and w = w(x, 2, t), whereupon Eq. (3), a linear partial 
differential equation, determines v = v(x, 2, ¢). 

The simplification achieved here is especially note- 
worthy because the problem represented by Eqs. (1), 
(2), and (4) and the boundary conditions for a given 
yawed cylinder is the same as for plane flow about the 
same cylinder at zero yaw. ‘Thus, the solutions for 
plane flow about any cylinder, can be carried over di- 
rectly to the yawed case, and the spanwise flow com- 
ponent v can be calculated, using such results, by inte- 
gration of a linear equation.* 

It is convenient to put Eqs. (1) to (4) in dimension- 
less form. According to the Prandtl argument re- 
garding orders of magnitude used above, the velocity 
components must have the formst 


u = 2, 
Vg(é, n, t’) (5) 
w = URe~ h(E, n, t’) 


where U and V are components of the free-stream veloc- 
ity vector (see Fig. 1) and ?’ = Ut/L, § = x/L, 7» = 


* It may be of interest to note that Eq. (3), for the spanwise 
velocity component v(x, z) in the steady-flow case, is identical 
with the equation for the temperature 7(x, z) in steady two- 
dimensional boundary-layer flow, provided that the Prandtl 
Number is unity and that viscous heating is neglected. [See 
Reference 2, p. 612, Eq. (35).] Thus, the distribution of span- 
wise velocity near a yawed cylinder is the:same as the distribution 
of temperature when its surface is maintained at a uniform 
temperature different from the stream temperature. This anal- 
ogy was called to the author’s attention by Prof. J. M. Wild, 
who also points out that the temperature distribution 7(x, z) 
is not affected by yawing the cylinder. 

Cf. Goldstein.? 


IIT. 


Re'"z/L, Re = UL/v, L being a characteristic linear 
dimension of the cylinder. Also, the pressure coef- 
ficient p’, defined as p/pU?, is determined by conditions 
at the outer edge of the boundary layer and is, in gen- 
eral, of the form 


= (6) 
The dimensionless form of Eqs. (1) to (4) is therefore 

fu + fhe + Why = + Sm (7) 

+ fee + = (8) 

pb,’ = 0 (9) 

fe +h, = 0 (10) 


where subscripts indicate partial differentiation. 


II. Discussion 


It is clear from the simplified form of the boundary- 
layer equations above that yaw of an infinite cylinder 
does not affect some of the important characteristics of 
its laminar boundary layer. In particular, the phe- 
nomenon of separation of the chordwise flow component 
from the surface, which occurs when f, vanishes at the 
surface (7 = 0), must appear at the same value of & 
as in plane flow. The result that the separation point is 
independent of Reynolds Number is known for two- 
dimensional laminar flows;? it can now be seen that it 
must be independent of yaw as well. 

It also seems possible to speculate on the effect of yaw 
on the phenomenon of transition. Apparently, transi- 
tion is ultimately produced by instability of the lam- 
inar layer.* Assuming that the yawed cylinder is 
subjected to the same flow disturbances, as might be 
the case in a homogeneous turbulent stream, one may 
conclude that the same instability must appear in the 
u and w velocity components of the yawed flow, at the 
same value of Re. Since Re is defined as UL/v, this 
means that the Reynolds Number based on the perpen- 
dicular velocity component controls the phenomenon 
of transition on a yawed infinite.cylinder. 


YAWED FLAT PLATE 


The simplest case of the laminar boundary layer is 
that of steady flow across a flat plate aligned with the 
flow direction, so that no pressure gradient appears. 
The two-dimensional problem was treated by Blasius,* 
who calculated the velocity profile that is usually known 
by his name. Using the equations above, one can now 
investigate the effect of yawing the leading edge of the 
plate, so that it is no longer perpendicular to the stream. 

The differential equations for this case are: 


+ = fan (11) 
Sg: + hg, = (12) 
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The solution of Eqs. (11) and (13) is the Blasius solu- 
tion : 
f=F(), h= — (14) 


where r = = U/ox. F'(r) is tabulated and 
plotted in several references.” * Now it is clear that, 
since the boundary conditions on f and g are the same 
—i.e., they both vanish for » = 0 and have the value 1 
at» = © for all §— and, since Egs. (11) and (12) are 
identical in form, the solution of Eq. (12) is: 


g=F(n) =f (15) 


Thus, « and v are everywhere proportional to one 
another, the constant of proportionality being U/V, the 
tangent of the angle of yaw of the leading edge. Hence, 
the conclusion is reached that the boundary-layer flow 
is not affected by yaw of the leading edge, and the 
boundary layer streams continually in the direction of 
the free stream.* 


IV. A CERTAIN CLAss OF YAWED CYLINDERS 


Prandtl has pointed out* that, in the steady two- 
diniensional case, an external-flow velocity of the form 


= ayx + (16) 


is sufficient to produce boundary-layer profiles of the 
“separation”’ type—i.e., having negative values of 
du/Oz at the surface ( = 0) — downstream of a certain 
location called the separation point. The parameters 
a, and a3 can be chosen to fit the velocity distributions 
of a variety of cylinders, at least in the regions upstream 
of their separation points. 

In accord with the assumption that p’ is a function 
of — only (in the steady case), it is clear that Eq. (16) 
is equivalent to 


um = UE — &?) (17) 


tified with the distance (—a,/a3)'/* between the loca- 
tions where “ vanishes and that U is a value propor- 
tional to, but not necessarily equal to, the speed of the 
undisturbed stream. Since the dimensionless surface 
velocity “;/U in potential flow is a given function of £ 
for any given cylinder, it appears that Eq. (17) is 
actually as general as Eq. (16). 

One can therefore adopt Eq. (17) to represent the 
external-flow conditions and pressure distribution for- 
ward of the separation point for a class of cylinders. 
The ratio of U to the normal stream speed varies from 
cylinder to cylinder. 

For this two-dimensional steady case, the boundary- 
layer velocities have been calculated.?»* The results 


* Note added in proof: In a new Technical Nolte of the N.A.- 
C.A. (No. 1402), R. T. Jones has considered the effect of yaw on 
the laminar boundary layer and has reached this same conclusion. 
He also points out the conclusion of Section II herein—i.e., that 
the flow in the xz plane is unaffected by yaw. 


provided only that the characteristic length L is iden-. 
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are given in the formt 


u = axfi'(n) + (18) 
w= —Vay filn) — 3asV — 

or 

f — + 

h = —fi(n) + 3&fs(n) — ... (19) 


The functions fi(), fi’(n), fa(n), etc., are tabulated in 
references 2 and 3. 

The separation point is given by fn(n = 0) = 0, or 
— = 0.69 Substituting Eqs. (19) in Eq. (8), one 
obtains the differential equation for g for the analogous 
case of the yawed cylinder: 


En — + lg: + 
[fi(n) — 38fs(n) + ....]g, = 0 (20) 


It is assumed that g will be found in the form 


g = 


n=) 


(21) 
Equating coefficients of powers of & one finds the fol- 
lowing set of ordinary differential equations: 


Go" + fiGo’ = 0 
Gy" + fiGi’ — fi'G: = 0 


G,” + fiGe’ 2f:'Ge 3fsGo’ (22) 
G3" + fiG;’ oF 3fi'Gs 3fsGr’ fs'Gi 
Gs" + fiGs’ — 4fi'Gs = 3fsGe2’ — 2fs’Go-— SfsGo’ 
etc. 
The boundary conditions are: 
G,(0) = 0, for all 
G,(~) = 0, forn 2 1 (23) 


G(~) = g(~) = 1 


Integration of the Differential Equations 


The appropriate solution of the first of Eqs. (22) is 
Go = ¢ fo” exp (— S'filn)dn)dn (24) 


This has been evaluated numerically, using the values 
of f; tabulated by Prandtl;* the value of the constant 
of integration c is found to be 0.570 The func- 
tion Go(n) is plotted in Fig. 2. 

It has been established that the differential equa- 
tions corresponding to odd values of nm in Eqs. (22) have 
no nontrivial solutions for the given boundary condi- 
tions. This is, of course, consistent with the condi- 
tions of symmetry of the present problem. 

The equations for n = 2, 4, 6, ... can be integrated 
numerically. The process is facilitated by transform- 
ing to a new independent variable which runs from 0 
to 1 when 7 runs from-0 to ©; for example, ¢ = 1 — 
e". The numerical integration can then be performed 


+ These results are not expected to apply in the neighborhood 
of the stagnation point, x = 0, where the boundary-layer ap- 
proximations are not valid. 
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Fic. 2. The functions Go(n), G2(n), Gi(n), and Ge(n) em- 
ployed in calculating the boundary-layer flow on a class of yawed 
cylinders. 
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Fic. 3. Profiles of the chordwise and spanwise velocity com- 
ponents at three locations in the boundary layer of a yawed 
cylinder. 


by subdividing the interval 0 S ¢ < 1 into N subinter- 
vals and approximating to the differential equation by 
a difference equation to be satisfied at N points. This 
procedure has been adopted here, with V = 50, for the 
numerical calculation of the functions Go, Gy, and Gs. 
The results are plotted against » in Fig. 2. 


Numerical Results 


Using the functions plotted in Fig. 2, one can calcu- 
late the spanwise flow component v with satisfactory 
accuracy within the region where the theory is ex- 
pected to hold—i.e., from a location somewhat aft of 
the leading edge to the location of separation of the 
chordwise flow, § = 0.69 .... (or perhaps to a small 
distance behind this location). In Fig. 3 are presented 
some typical profiles of the velocity components u and 
v; these have been combined in Fig. 4 to show the 
boundary-layer flow in both direction and magnitude. 

In Fig. 5 the “‘limiting streamline’ at z = 0 has been 
drawn in comparison with the corresponding potential- 
flow streamline of the outside flow. The “limiting 
streamline’ is the curve whose direction everywhere 
coincides with that of the (vanishing) fluid velocity at 
the surface; thus its equation is 


dy/dx = lim (v/u) 


The scales chosen to plot the streamlines in Fig. 5 are 


v =4 
0 os 10 os 


N \ 
| 
4 
B= 0.5 —=0.7 

Fic. 4. ‘Vector diagrams showing the direction and magnitude 

of velocity at three locations in the boundary layer of a yawed 


cylinder. The parameter 7 defines the dimension from the sur- 
face of the cylinder. 
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Fic. 5. The streamline of the potential flow and the limiting 
streamline at the surface, for flow over a yawed cylinder. 


arbitrary and have simply been selected to show the 
streamlines clearly, since neither the ratio V/U nor the 
ratio of U to the free-stream normal component has 


‘been specified. 


Other streamlines of the boundary-layer flow are more 
difficult to determine, since they require evaluation of 
w. It appears from available numerical data that, in 
plan view, they will lie between the two extreme cases 
drawn in Fig. 5. 
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Optimum Number of Webs Required for a 
Multicell Box Under Bending’ 


GEORGE GERARDt 


Republic Aviation Corporation 


SUMMARY 


In multicell construction that may be used on the wing and tail 
surfaces of high-speed aircraft, the relatively thick compression 
surface is stabilized by a series of webs. It is necessary to deter- 
mine the correct combination of number of webs and compres- 
sion cover thickness which will result in the design of minimum 
weight. 

Theoretical results are obtained in which the optimum number 
of webs is determined to be a function of only the structural thick- 
ness ratio. Nondimensional design charts are presented from 
which optimum conditions can be obtained. In addition to op- 
timum conditions, an investigation was made to determine the 
weight penalty involved in using a number of webs other than 
optimum. It was found that for a slight weight penalty it is 
possible to reduce the number of webs required and still maintain 
an efficient design. 


INTRODUCTION 


| he CERTAIN TYPES of high-speed aircraft wing and 
tail surfaces, it may be desirable to utilize multicell 
construction in which relatively thick skins are stabil- 
ized by a series of webs. For the multicell box shown 
in Fig. 1, acting under a pure bending moment, the 
principal role of the webs is to displace the effective 
bending material as far from the neutral axis as possible 
and to divide the compression cover into a series of 
panels of small width-to-thickness ratio to attain high 


‘critical compressive stresses. 


The first role of the webs can be simply met by only 
one extremely thin web along the vertical centerline of 
the box, forming, in effect, an extremely wide flange 
I-section.t The second role of the web, that of sup- 
porting the compression cover against instability, re- 
quires a series of webs to divide the cover into panels of 
high critical stress. This can be achieved with various 
numbers of webs, and, consequently, it is necessary to 
establish web and compression cover dimensions that 
will result in minimum structural weight for the load 
to be carried at a given spanwise section. . 


Received June 24, 1947. 

* Shortly after this investigation was completed in May, 1947, 
as part of a larger research program on thick skin structures, an 
N.A.C.A. report on the same subject was released. For complete- 
ness, this report is cited: Schuette, E. H., and McCulloch, J. C., 
Charts for the Minimum Weight Design of Multicell Wings in 
Bending, N.A.C.A. T.N. No, 1323, June, 1947. 

¢ Principal Research Engineer, in Charge of Structural Re- 
search Group. Now, Instructor, Guggenheim School of Aero- 
nautics, College of Engineering, New York University. 

t Note that the box is assumed to act under pure bending 


moment. 


It has been found from tests on heavy skin multicell 
boxes that instability of the thick cover is tantamount 
to failure of the box. The large loads introduced into 
the relatively thin webs at the onset of instability of 
the heavy compression cover are sufficient to cause com- 
plete collapse of the structure. Therefore, it is not 
reasonable to design the cover to buckle at limit load 
factor and expect the structure to retain structural in- 
tegrity until ultimate load factor, an additional load 
carrying capacity of 50 per cent. 

In addition, it has been determined! that the com- 
pression cover may be efficiently designed when the 
critical stress of the skin is equal to the ultimate applied 
load on the surface.’ Hence, under these considera- 
tions, efficient design of a multicell structure probably 
can be achieved when bending instability of the webs is 
attained simultaneously with compressive instability 
of the cover at the ultimate applied load.** 

It is conceivable that a design may be obtained 
utilizing two webs and a thick cover of relatively low 
critical stress or, in the other limit, a thin cover sup- 
ported by many webs which would result in a high 
critical compressive stress. From a weight standpoint, 
such a design would be prohibitive, since in the former 
case the cover would be unduly heavy while in the lat- 
ter the multiplicity of webs would result in excessive 


** In this analysis, the secondary compressive forces intro- 
duced into the webs by curvature of the multicell box under 
bending are neglected. 
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weight. It is evident that some arrangement inter- 
mediate to these two limiting conditions will result in 
a structure of minimum weight depending upon the 
number of webs chosen. 

Hence, the objective of this investigation is to deter- 
mine the optimum number of webs from a weight stand- 
point and to determine the parameters that influence 
this choice. In addition, since a minimum number of 
webs is degjrable from a fabrication standpoint, the 
weight penalty involved in departing from the optimum 
conditions will be investigated. 

Although the multicell box is considered to be under 
a pure bending moment for simplicity of analysis, it is 
felt that results obtained are significant when applied 
to an actual cantilever wing. Consideration of shear 
in the webs will probably result in a fixed increase in 
web material in all cases of multicell construction if the 
stresses are below the proportional limit of the material. 
The smallest increase in web material will probably ac- 
crue to the largest number of webs. In any case, the 
results of the investigation will indicate the range of 
optimum design for use in preliminary design work. 


SYMBOLS 


Refer to Fig. 1. 

= thickness, in. 

= width, in. 

modulus of elasticity, lbs. per sq.in. 
stress, lbs. per sq.in. 

Poisson’s ratio 

number of bays (x + 1 = no. of webs) 
solidity 

loading per chordwise in., Ibs. per in. 
instability coefficient, also K = kw?/12(1 — v?) 
solidity factor 


t 

b 
E 
v 
n 
z= 
Nz 
k 
? 


Subscripts 


cr. = critical 

$s = compression cover 
t tension cover 

w web 


ANALYSIS 


Assumptions 


For purposes of analysis, the wing section shown in 
Fig. 1 is idealized to the section shown in the lower part 
of that figure. In addition, the following assumptions 
are made: (1) the webs act under pure bending only; 
(2) the contribution of the webs is to support the cover 
only and take no part of the load; (3) the webs are 
essentially line elements with no caps at attachment 
to the cover; (4) the elements subject to instability are 
conservatively assumed to be simply supported; and 
(5) the box section possesses no chordwise chamber. 


Critical Stress Relations 


The critical compressive stress of the cover is taken 
to be 


Ger, = [ksr?E/12(1 — v?)](t,/b,)? (1) 


where k, has a value of 4 for the assumed conditions of 
simple support, large aspect ratio, and negligible cover 
taper. 

The critical stress in bending of the web is likewise 
taken to be 


Cory = [Rum?E/12(1 |(tw/bw)? (2) 


where k,, has a value of 25 when the web is simply sup- 
ported and of large aspect ratio. 

As discussed in the introduction, efficient design of 
the multicell box assumes the web critical stress to be 
equal to the cover critical stress. Consequently, by 
equating Eqs. (1) and (2) to meet the condition o,,., = 
Ccr, the result is obtained that the relative cover and 
web dimensions are 


(5/2) (tw/bw) 


or 


ty/ts = (2/5)(b./bs) (3) 
Solidity 

In calculating the weight of the elements of the multi- 
cell box section shown in Fig. 1 which are subject to 
inst&bility, it is convenient to utilize the concept of 
“solidity.’’ Solidity is defined as the ratio of the area 
of the structural elements comprising the cross section 
to the total area of the cross section enclosed by the 
outer contours. It is evident that the solidity bears a 
direct proportionality to the weight of the structural 
material. 

In considering optimum number of webs, the solidity 
of the idealized box section shown in Fig. 1 need include 
only the compression cover and the webs. The tension 
surface is designed from considerations of the ultimate 
tensile strength of the material and, consequently, for 
a given loading, N,, the thickness of the tension skin 
is constant regardless of the number of webs employed. 

Thus, consideration of the tension surface adds a con- 
stant value to the solidity of the structure and does not 
influence the optimum number of webs. The total 
solidity, =, of the box structure can be expressed as 


or 
nb, + + 1)Dutw te 
— 4 
3; + (4) 


The last term is the contribution of the tension surface 
to the solidity 


De 


and, being a constant, drops out in the development 
that follows. Therefore, the solidity of the elements 


subject to instability in Eq. (4) reduces to 
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TABLE 1 
Solidity Factor as a Function of Structural Thickness Ratio and Number of Bays 
n=1 n=2 n=3 n=4 n=5 n=6 n=7 n=8 n = 10 n = 12 
bw/nb, ¢ ¢ ¢ ¢ ¢ ¢ 
0.05 20.04 12.48 9.74 8.09 7.04 6.33 §.77 5.36 4.78 4.41 
0.10 10.08 6.35 5.04 4.28 3.83 3.54 3.34 3.22 3.11 3.10 
0.15 6.79 4.35 3.55 3.14 2.89 2.78 2.76 2.77 2.86 3.06 
0.20 5.16 3.40 2.87 2.62 2.53 2.54 2.59 2.69 2.98 3.38 
0.25 4.20 2.85 2.50 2.38 2.39 2.48 2.62 2.80 3.23 3.74 
0.30 3.59 2.52 2.29 2.27 2.37 2.54 2.74 2.99 3.56 4.20 
0.35 3.14 2.29 2.18 2.24 2.42 2.68 2.92 3.24 3.93 4.71 
0.40 2.82 2.14 2.12 2.26 2.50 2.79 3.13 3.51 4.33 5.24 
0.45 2.58 2.05 2.11 2.31 2.61 2.96 3.36 3.80 4.75 5.78 
0.50 2.40 1.98 2.12 2.38 2.74 3.15 3.60 4.10 5.17 6.34 
> t, | 1f/mb, 4 n+ 1/ty (5) To obtain the optimum number of webs which will 
wb Ln\ by n ts ' yield minimum solidity, it is only necessary to differen- 


From Eq. (3), however, a relationship for ¢,,/t, has been 
found. Substituting this into Eq. (5), 


ts 2(n + 1) 
Loading 
In Eq. (6), the structural thickness ratio 5,,/nb, is 
known and n, the number of bays, must be an integer. 
Only ¢, is not known but this term bears a simple rela- 
tionship to the applied loading on the section. The 
loading, N,, is the load on the cover per chordwise inch 
and, according to reference 1, 


N, = (7) 


Fer.sts 

Since the expression for ¢,,, is given by Eq. (1), 
N, = K,Et,*/b,? 

or 


t, = WN,b2/K;E (8) 


Optimum Conditions 


By combining Eqs. (8) and (6), the following expres- 
sion for solidity is obtained in terms of the loading, 
N,, structural chord, nb,, structural thickness ratio, 
b,,/nb,, and number of bays, n. 


The cube-root term of Eq. (9) is merely a multiplying 
factor for determining the solidity since nb, is the struc- 
tural chord. By transposing this term, it is possible to 
obtain a “‘solidity factor,’”” ¢, in terms of structural 
thickness ratio and number of webs only. 


N, 
K,Enb, 


Nz 


¢ 


—1/; 
K 


tiate Eq. (9) and set it equal to zero, O2,,./0n = 0. 
Performing the differentiation leads to the following 
expression for the optimum number of bays: 


n(4n + 1) = 5 (nb,/by)* (11) 


It is highly significant that the optimum condition 
depends only upon the structural thickness ratio as 
indicated by Eq. (11) and is independent of the critical 
compressive stress of the cover. 


DISCUSSION OF RESULTS 


Based on the assumptions of the analysis, the follow- 
ing analytical results have been obtained. The solidity 
factor, ¢g,.is a function of the structural thickness ratio 
and number of bays only and is independent of the 
loading and structural chord. Apparently, it depends 
solely upon the multicell geometry. The expression 
for solidity factor is 


In Table 1, values of solidity factor have been com- 
puted for various thickness ratios and numbers of bays. 
These data are plotted on Fig. 2 and show the number of 
bays required for minimum solidity as a function of 
structural thickness ratio. The lower envelope of this 
curve follows the equation for the optimum number of 
bays. 


(11) 


In current designs of most two-spar or multispar 
wing or tail surfaces, the structural chord is between 
0.5 to 0.6 of the unswept aerodynamic chord. Hence, 
a simple method of using Fig. 2 is to double the aerody- 
namic thickness ratio to estimate the structural thick- 
ness ratio. Knowing the latter, the optimum number 
of bays and solidity factor can be found directly from 
aerodynamic data. 

While Fig. 2 gives the optimum number of bays based 
on the assumptions of the analysis for minimum weight, 
fabrication difficulties may preclude use of optimum 
conditions in design. Consequently, it is necessary to 


n(4n + 1)=5 (nb;/b.)* 
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Fic. 2. Solidity factor as a function of structural thickness ratio. 


consider the weight penalty involved in using less webs 
than is indicated by the optimum. 

Since weight is directly proportional to the solidity 
factor, a certain percentage increase in solidity factor 
will result in a corresponding increase in weight. On 
Fig. 3, the data of Fig. 2 are cross-plotted to indicate 
the weight penalty involved in using a number of webs 
other than optimum. For a relatively small weight 
increase of 5 per cent, Fig. 3 indicates the results shown 
in Table 2. 


TABLE 2 
Number of 
Structurai Number of Bays* for 
Thickness Bays* for 5 Per Cent In- 
Ratio Optimum crease in Weight 
0.10 11 8 
0.15 7 5 
0.20 5 4 
0.30 4 3 


* To the nearest integer. 


Fig. 3 is particularly useful when the structural thick- 
ness ratio may vary spanwise in a given wing with the 
result that a compromise must be made in selecting the 
best design. 


CONCLUSIONS 


Based on the assumptions of the analysis, the follow- 
ing results have been obtained: 

(1) A theoretical expression for the optimum num- 
ber of webs required in a multicell box for minimum 
weight has been derived in which optimum conditions 
depend only upon the structural thickness ratio of the 


2 4 


Fic.3. Weight penalty for number of bays other than optimum. 


section. Use of Fig. 2 permits determination of the 
optimum number of webs. 

(2) The weight penalty involved in using a number 
of webs other than optimum has been determined from 
the theoretical expression. For a slight increase in 
weight, it is possible to reduce the number of webs re- 
quired and still maintain an efficient design. 
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On the Kinematics of Turbulence 


F. N. FRENKIEL* 
Cornell University 


SUMMARY 


A short review of different correlation coefficients and spectra 
is presented. Some new relations between these characteristics 
are given for homogeneous and isotropic turbulence. The ‘‘form 
of a turbulent particle’ has been given from the point of view of 
correlation. 

The study of numerous experimental correlation curves shows 
that, where their second moments are represented as functions of 
their first moments, the resulting points lie near a straight line. 

The correlation curves may be represented by functions of the 
genera] form 


R(r) = exp(—k|r|)®(r) 
or 
exp(— kr?) ®(r) 


Several examples are given of correlation functions fitted to ex- 
perimental results. The use of such functions simplifies the com- 
putation of spectra or of the transverse and longitudinal correla- 
tion curves, one from the other. This may be useful in studying 
the history of spectra of turbulence behind grids and in the 
boundary layer. 


INTRODUCTION 


— STATISTICAL THEORY OF TURBULENCE has 
been the object of many investigations and 
publications. The papers of Th. von Karman, H. L. 
Dryden, and G. I. Taylor in this JouRNAL and else- 
where have presented various aspects of the subject. 
In a more recent paper, Dryden! has made a complete 
review of the theory. 

The present paper is concerned with the kinematics 
of turbulence and, more particularly, with the repre- 
sentation of correlation curves and spectra of turbu- 
lence by empirical functions. 

It seems to the author that it may be useful to make 
a clear distinction between the different correlation 
coefficients by employing appropriate notations. Simi- 
lar to the spectrum of turbulence at a fixed point (G. I. 
Taylor) and the spectrum studied in the Lagrangian 
system (J. Kampé de Feriet), the longitudinal and 
transversal spectra are used here. 

The relations between the different characteristics 
of isotropic turbulence are completed by some new 
equations. 

The results of the experiments of H. L. Dryden and 
his collaborators, L. F. G. Simmons, A. A. Hall, and 
D. C. Macphail, are here represented by empirical 
functions. When either the transversal or longitudinal 


Received May 26, 1947. 
* Previously at the Groupement Francais pour le Développe- 


ment des Recherches Aéronautiques. 


correlation function is given, it is easy to compute the 
other for a flow of isotropic turbulence. It may then 
be verified whether the experimental correlation points 
correspond to the computed curve and ascertained 
whether the turbulence seems to be isotropic. The com- 
putations of the spectrum from the correlation curve, 
or the converse may also be made, employing the repre- 
sentative functions. A more detailed report on the 
empirical correlations and spectra functions is included 
in another paper. 


TURBULENT AND MEAN VELOCITIES 


In a turbulent fluid flow the instantaneous velocities 
vary in time and in space in a manner that is, or seems 
to be, disorderly. The only means of direct characteri- 
zation of the turbulence is then by use of statistical 
values. 

Let V(x, y, 2, t) be the instantaneous vectorial velocity 
at a given point P(x, y, z) at a given instant t. The 
quantity 


m(V) = = lim (1/2T) V(x, 2, s)ds 


in which the integral is taken in a vectorial manner, is 
called vectorial average velocity of the fluid at P. 
The turbulent vectorial velocity is given by the rela- 
tion 
V'(x, y, 2, t) = V(x, y, 2, t) — V(x, y, 2) 


The instantaneous vectorial velocity V may be con- 
sidered as the vectorial average velocity of the mole- 
cules that .are present in a space-time domain of an 
elementary volume v, enclosing the point P and in an 
elementary time interval 7, containing the instant ¢. 
The domain 7,v; is assumed to be sufficiently large with 
respect to the molecular scale so that this average may 
have a definite value, and small enough with 
respect to the fluid motion, to allow the neglect of the 
turbulence in 

Let Oxyz be orthogonal coordinate axes such that Ox 
is parallel to the direction of the mean velocity V. De- 
note by u, v, w the projections of V in the three direc- 
tions. The scalar mean velocity components are given 


by 
1 t+T 
M,(u) = = lim u(x, y, 2, s)ds 


=d=0, 
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TABLE 1 
Study of turbulence is made with the: = | Scale of 
Turbulence 
Components of Correlation Coefficient (Length or Spectrum of Turbulence 
the turbulent At two points | Measured at Time of — 
velocity placed on the instants * Notation Name | Correlation) | Notation Name 
Parallel to a di- A line paral- Simultane- (r) Axial corre- Lyk (w) 
rection A lel to the di- ously at lation 
rection A an in- 
Perpendicular to and at a stant (r) Tangential Lat (@) 
a direction A distance correla- 
equal tor tion 
Ra(r) = Ra"(r) Correlation La fa(w) Spectrum of tur- 
along A bulence along A 
A line parallel 
to Ox and Rz(x) = Rz (x) Longitudinal Lz S2(w) Longitudinal spec- 
at a_ dis- correlation trum of turbu- 
tance x lence 
Parallel to the, A line parallel Ry(y) = Ry*(y) Transverse Ly Transverse spec- 
direction of to Oy and correlation trum of turbu- 
the mean ve- at a_ dis- lence 
locity chosen tance y 
as axis Ox 
At point Ri(h) = Re“ (h) Correlation Li St(w) Taylor spectrum 
fixed in re- at a fixed . 
gard to At two in- point 
Oxys stants 
andi +h 
At a_ point Rm(h) = Rm*(A) Correlation Lan Spectrum follow- 
moving with mean ing mean mo- 
with mean motion tion 
velocity 
with respect 
to Oxyz 
Following the Rit(h) = Ret“ (A) Diffusion Kampé de Feriet 
same turbu- correlation spectrum 
lent particle 
Parallel to axis Rar(r), Rz(r), La”, Le", Ly”, fa*(w), fe"(@), 
Oy etc. ete. etc. 


To the vector of turbulent velocity V’ correspond 
three turbulent components, u’ = u — a, v’ = v, and 
w’ = w, and asa result of the definition of the turbulent 
components, their averages will be equal to zero: 


=wy’ = 0. 


CHARACTERIZATION OF TURBULENCE 


To characterize the turbulence of a fluid flow, the 
expressions most commonly used are: , 

(1) Mean-square values of the turbulent velocity 
components: The ratio of the root-mean- 
square value of the turbulent velocity to the mean 
velocity # is called the intensity of turbulence. Three 
components of this characteristic may be defined: the 
intensity of longitudinal turbulence, T, = Vv. u!2/%, 
and two intensities of transversal turbulence, 


T, = and T, = 


_ (2) Correlation coefficients, such as 


u’ (x1, V1, 21, t1)U" (Xe, Yo, te) 
Re 
(x1, My 21, th)? [u! (22, 2, Za, te)]? 


between parallel components of the turbulent velocities 
at two points of the domain space-time, or correlation 
coefficients between two perpendicular components, 
such as u’, v'; u’, w’ or v’, w’, at the same point or at 
two points of the space-time domain. The correlation 
may be studied in the Eulerian system, using the veloci- 
ties at fixed points with reference to the coordinate 
axes with no identification of the actual turbulent par- 


ticles passing these points. In relation to turbulent dif- 
fusion it is also necessary to use the correlation coef- 
ficient of the Lagrangian system. In this case the veloc- 


ity fluctuations are studied by following the turbulent — 


particles along their paths. It is sometimes useful to 
study the correlation following the mean motion of the 
fluid and thus to characterize correlation coefficients 
in a pseudo-Eulerian system. 


(3) Scales of turbulence are obtained by integrating 
the correlation curve with reference to one of the param- 
eters x, y,zorh. The scale of turbulence, which may 
be called correlation length or correlation time, accord- 
ing to the parameter that is studied, is equal to L = 
So” R(s)ds. 

(4) Spectra of turbulence f(w), represent the turbu- 
lent energy distribution as a function of the cyclic fre- 
quency w = 2rn (n = frequency). f(w)dw is defined as 
the amount of turbulent energy for which are responsi- 
ble the frequencies between w and w + dw. These 
spectra may be obtaizied by harmonic analysis of the 
fluctuations of u’(x, y, 2, f) given as a function of one of 
the parameters. 


Some of the characteristics of the turbulence are 
given in Table 1, and Fig. 1 presents schematically the 
correlation coefficients as functions of the distance. 
The notations employed may not seem simple, but it is 
believed that it will be desirable to use different nota- 
tions for the different characteristics of turbulence and, 
thus, to avoid confusion between the correlation coef- 
ficients and between the spectra used in the statistical 
study of turbulence. 


Th 


will 


4 
> 
i 
the ti 
2 
the c 
i 
lim ( 
i then 
Euler 
4 
ii equal 
4 As 
comp 
spect 
from 
ll 
| 
; 
G. 
y 
ment 
Nati 
Sti 
Kam 


trum 


‘ollow- 
1 mo- 


Feriet 


THE KINEMATICS OF TURBULENCE 


59 


Fie. 1. 


The turbulence is called homogeneous in y when 


(1) 


i.e.—when a translation of the axes does not change the 
average values of the fluctuations. When, in addition, 
the time averages are equal to the space averages, as in 
the case when 


lim (1/2T) fits" W(x, y, 2, s)ds = lim (1 /2X) 


¥(x, y, 2, t) = F(x, y, 2) = const. 


y, 2, Dds = lim (1/2V) W(x, s, t)ds 
Yoo ‘ 


(2) 


then the space correlation coefficients in the 
Eulerian system and in the Lagrangian system are 
equal.* 

Assuming that the turbulent velocities are small in 
comparison with the mean velocity, G. I. Taylor‘ found 
that the longitudinal correlation coefficient and the 
spectrum of turbulence at a fixed point can be calculated 
from one another by Fourier’s transformation: 


fle) = (2/nit) va (“) R.(s)ds 


G. I. Taylor has verified these equations for experi- 
ments made by L. F. G. Simmons and C. Salter in the 
National Physical Laboratory. 

Studying turbulence in the Lagrangian system, J. 
Kampé de Feriet® found similar relations 


(3) 


fu(o) = “cos (ws)R(s)ds 


will be 


Correlation coefficients of the second order between simultaneous velocities at two points. 


Riz(h) = Jo” cos (sh)fi(s)ds (4) 


without the assumption of u’ < @. 

By using, in an analogous way, the definition of the 
longitudinal spectrum of turbulence given above, it can 
be found that 


f(w) = 2 cos (=) R.(s)ds 


R,(x) = cos (=) fr(s)ds 


without the restriction concerning the magnitude of 
the turbulent velocities. . 

In the case for which G. I. Taylor has verified the 
equations, it is then possible to assume that 


filw) = 


(5) 


and 


Ri(h) = (6) 


This equality does not necessarily exist in every tur- 
bulent flow. 


IsoTROPIC TURBULENCE 


In a flow of homogeneous and isotropic turbulence, 
the averages of the turbulent velocities are independent 
of a translation or of a rotation of the axes. In that 
flow = = wand = y’w’ = = 0. von 
Karman® has demonstrated that between the two 
coefficients of correlation there exists the relation 


RA(r) = RX(r) + (1/2)r[dRX(r)/dr] (7) 
To this equation there corresponds the inverse relation 


Ri(r) = (2/r*) fo" sRa(s)ds (8) 


In particular, for the correlations between the components parallel to the mean-velocity direction, the relations 
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By integrating the two sides of the first of these two relations, noting that the condition lim n fe(w) = 0 leads to 


fin 3 x -R 2(x) = 0, it will be found that the longitudinal scale of turbulence is equal to twice the <a scale of 


L, = 2L, (10) 


The correlation coefficient along A may be given as function of the longitudinal and transverse correlation coef- 
ficients by 


= = Re) Rl) (uy 


where r = V x2 + y?, If a is the angle between the direction A and the x-axis, then 
R,(r) = R,(r) cos? a + R,(r) sin? a 
After integrating the last equation from zero to infinity and applying Eq. (10), it will be found that 
Ly = L,(1 + cos? a) (12) 


Fig. 2 represents the curve that corresponds to Eq. (12). If one wishes to imagine the “turbulent particle’’ 
from the point of view of the turbulent components parallel to the direction Ox, then the figure formed by revolving 
this curve about 0x may be used as a convenient picture. 

To the two Eqs. (9), which give the relation between the longitudinal and the transverse coefficients of correla- 
tion, there correspond the two following relations between the longitudinal and transverse spectra of turbulence: 


as a function of C = Vo, starting from the relation 
= —(d?R,(0)/dx?) 
and applying the equation 
(1/u?) fo” s*f.(s)ds = —(d°R,(0)/dx*) 
It will be found that 
C= (15) 


(11a) 


and f;(w) 


An analogy is often made between the spectrum of 
turbulence and the spectrum of light. This comparison 
may be continued by defining a factor that presents 
some analogy with color. It is proposed to designate 
as “color of turbulence” the standard deviation of the 
spectrum of turbulence. The square of the “color of 
turbulence’ will be given by 


= w= s'f.(s)ds (14) 


The quantity A, introduced by G. I. Taylor and repre- 
senting the dimension of the smallest eddies, responsible 
for the dissipation of turbulent energy, is determined 


REPRESENTATION OF CORRELATION BY FUNCTION 


The possibility of representing the correlation curve 
or the spectrum of turbulence by a function will sim- 
plify the study of turbulence. It is evidently possible 
to use general methods of interpolation and to give a 
representation of the experimental results by an alge- 
braic or trigonometric polynomial. Nevertheless, this 
method is not wholly advisable. It seems important 
to reckon with the very particular form of the experi- 
mental curves. To represent the correlation curves, 
Dryden’ used a function of the form 


R,(y) = exp(—|y|/L,) 


A similar function was found theoretically by Doob*® 
to represent the correlation in Brownian movement. 
On the other hand, the experimental curves bring to 
mind the normal error function of Gauss, and it seems 
of interest to study its application as well. 


Fic. 2. Form of a “turbulent particle’ considered from the 
point of view of the turbulent es parallel to the direction 
of the mean velocity. 
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In Fig. 3 the abscissa and the ordinate represent the 
first and second moments of certain experimental cor- 
relation curves published by different authors.* ” -" 
for the computation of these moments, only the parts 
of the curve for which R is positive have been taken 
jnto account—i.e., the values are given of 


1 Po 
x 


and 
= 
x 
where 
x = A” Ga(s)ds 
Ra(p) = Ra(r/Ly) = Ra(r) 
po = ro/Ly 


and 7) represents the abscissa r of the experimental cor- 
relation curve at the point where R, = 0. It is noted 
that the representative points of these moments are 
grouped near a straight line. The points corresponding 
to the two functions R,(p) = exp(—|p!) and R,4(p) = 
exp(—p?/4) are near the same line. ' 

This suggests the use of functions having the form 


R(r) = exp(—ér”™) 


with 1 < m & 2 for the representation of the correlation 
coefficients. As the use of such functions has not 
proved easy, it has been preferable to employ functions 
of the form exp(—k|r|)®(r) and exp(—kr*)®(r), in 
which #(7) is an algebraic or trigonometric polynomial. 

Four conditions, necessary but not sufficient, are put 
down for the correlation functions; viz., 


-1¢Rr)< +1, SL” R(s)ds > 0 
lim R(r) = 1, lim R(r) = 0 


(16) 


The verification of these conditions limits the extent of 
values that may be given to the arbitrary coefficients 
of ®(r). In another publication? are given the limits 
corresponding for a certain number of functions. 


COMPARISON WITH EXPERIMENTAL RESULTS 


Some examples of experimental results are compared 
below with chosen functions for correlation coefficients 
or for spectra of turbulence. Figs. 4 and 5 represent 
experimental results obtained by Dryden and his col- 
laborators’ at the National Bureau of Standards, for the 
transverse correlation coefficient measured at a distance 
of 40 mesh lengths behind two grids of 5-in. and 3'/4- 


in. mesh sizes. The curves drawn on both figures cor- 


respond to the equation 
®,(n) = [Ao + Ai cos (men)] exp(—cln)) 
where 7 = y/L, and &,(n) = R,(y). There are two 
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Fic. 3. The relation between first and second moments of ex- 
perimental correlation curves. 
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Fic. 4. Transverse correlation behind a grid of 5 in.” 
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Fic. 5. Transverse correlation behind a grid of 3/, in.” 
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Fic. 6. Transverse and longitudinal correlations at 85 mesh 
lengths behind a grid of !/, in.® 


B=0,25 
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Fic. 7. Transverse and longitudinal correlations at 27.5 mesh 
lengths behind a grid of 3 in.‘ 


arbitrary coefficients, m and Ap, in this equation. A, c 
and the limits between which m and A» may vary are 
found by applying the conditions given in Eqs. (16). 
For the actual function the result will be 


A, = 1— Ag, c = Aym? + 1/m? +1 
and 


—1/m < Ao < 
[(5m? + 8) 16m* + 80m? + 64]/9m? 


The result of the experience with the 5-in. mesh grid is 
represented by the function with coefficients m = 0.5 
and Ao = 2.75. For the tests behind the 3!/,-in. grid, 
the values are m = 0.5 and Ap = 2. 

For some experiments, satisfactory results have been 
obtained with functions of the general form 


= [1 + 
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Fig. 6 shows, as a function of » = r/L, (n is either «/L, 
or y/L, as the case may require), the experimental] 
points obtained by Hall’ for the correlation at 85 mesh 
lengths behind a grid of !/;-in. mesh. The round points 
represent results of the measure of the transverse cor- 
relation coefficient ®,(y) and the square points the 
longitudinal correlation coefficient ®,(n). The @, 
points are fitted by a function 


= [1 + 


with A; = 0.95 and ¢ =I + A;. When the second of 
Eqs. (9) is applied to compute a corresponding &,- 
function in a flow of homogeneous and isotropic turbu- 
lence, then it is found that 


1 


[(1 + 2A1)(1 + + 


with » = 2& The curve corresponding to this function 
well represents the experimental points of the longi- 
tudinal correlation coefficients. The turbulence of the 
flow in which Hall made the measurements could there- 
fore have been nearly homogeneous and isotropic. 

The results, which have already been given in Fig. 4, 
may be also represented by the functions 


Ry(n) = [1 + Are|n| + 


with A; = 0.10, As = 0.25, andc = 1 + A; + 242. 
Some results may be represented by the relations of 
the general form 


R(p) = A, exp(—c, |p!) 


Typical are the experimental correlations measured by 
L. F. G. Simmons at a distance of 27.5 mesh lengths 
behind a grid of 3-in. mesh. These experiments, in 
which the longitudinal and the transverse correlation 
coefficients were determined, have been already studied 
by Taylor.‘ In Fig. 7, the points giving ®,(£) are rep- 
resented by a curve corresponding to the function 


= A exp(—clé|) + B 


where A = 0.25, 6 = 0.25,c = A + B/8B, and B = 
1-A. 

When the transverse correlation function is com- 
puted, it is found 


= — (1/4)e\n|] exp[—(1/2)eln]] + 
Bll (/4)Be|n|] exp[—(*/2)8c|n| 


The corresponding curve is drawn on the Fig. 7 as a 
function of — = 7/2. It is seen that the two curves 
and R,(é) represent the experimental points. 
The turbulence of the flow in which the experiments 
were made can certainly have been homogeneous and 
isotropic. 

Fig. 8 represents the results of experiments made by 
Macphail.!? The square points represent the correla- 
tion coefficient ®,’ between the components of the 
turbulent velocities perpendicular to the direction of 
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the mean motion. This correlation is fitted by the 


curve 
= A exp(—eln’|) + B exp(—Beln’|) 

where 7’ = r/L,’ and with the coefficient 6 = 0.50 

and A = —0.9. 


The round points represent the correlation coef- 
ficients ®.° and ®,, which for isotropic turbulence 
would be equal to ®,”. The curve ®,” = &.’ = ®, has 
been computed with an equation similar to Eq. (7). 
The positions of the experimental points relative to this 
curve seem to show that the turbulence of the flow in 
which these experiments were made was not homogene- 
ous and isotropic. 

In a similar manner it is sii to represent a spec- 
trum of turbulence by a function. The representative 
equation may be given directly or as a Fourier trans- 


' form of a chosen correlation function. Fig. 9 gives the 


points of a spectrum f,(w) measured by L. F. G. Sim- 
mons and C. Salter and studied by Dryden."* The 
curve drawn on this figure corresponds to the correla- 
tion function 


Rr) = A exp(—c?r?) + B exp(—’c?r?) 
in which r = h/L, A = 0.6, 8 = 0.15, and where 
B=1 —Aand 


= (Vx/28)[A(6 — 1) + 1] 


are obtained as a result of the conditions (16). 

The corresponding spectrum’ function, calculated 
with an equation similar to Eq. (4), has the same gen- 
eral form as the correlation function: e 


1 2? B 2? 


with 2, = wl, and = fi(w)/Li. 


CONCLUSIONS 


The use of representative functions to characterize 
correlation is particularly useful to facilitate the calcu- 
lation of Fourier transformations [Eq. (3)] and to 
carry out the numerical calculations required to verify 
the isotropy of turbulence. Graphical methods for 
these computations are not precisé and in general are 
laborious. For the experiments of Hall (Fig. 6) and 
especially for those of Simmons (Fig. 7), von Karman’s 
equation is verified. On the contrary, the results of 
Macphail (Fig. 8) seem to have been obtained in a flow 
of anisotropic turbulence. . 

The representative functions must be carefully 
handled to avoid errors, The experimental results 
usually cover time or distance intervals that are not 
very large. The difference between the correlation 
coefficients resulting from the choice of the function 
and the real coefficient may sometimes be important 
in intervals for which the results are not available. As 
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Fic. 8. Correlation according to Macphail’s experiments at 30 
- mesh lengths behind a grid of 2 in.” 
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the use of the function implies a preliminary integra- 
tion, these differences may add up and give an impor- 
tant error. The use of the experimental value of L,/L, 
as a criterion of isotropy must then be employed care- 
fully. The conclusion given by Macphail, as a result 
of his experiments, is quite contrary to that given here 
above, as far as isotropy is concerned. It may be no- 


ticed that the evaluation of L,/L, by graphical extra- 
polation seems to confirm the nonisotropy of these ex- 
periments. It is possible—but does not seem prob- 
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able—that the isotropy of Macphail’s turbulence 


“could be confirmed by the methods used here, if 


more complicated representative functions were em- 
ployed. 


It is believed that the use of representative functions 
may help the study of the transformation of the spec- 
trum of turbulence in the flow. A turbulence spectrol- 
ogy will certainly be more difficult than the spectrology 
of light. A filter traversed by the light passes a charac- 
teristic band of frequencies that are conserved in the 
spectrum. If the turbulence behind a grid is studied, 
the large eddies break up into small ones which, con- 
tinuing the process, finally dissipate in molecular agita- 
tion. The initial frequencies generated by the grid are 
not conserved as in the case of light. Nevertheless, 
one may hope to obtain some information about the 
prior history of the flow from the study of a turbulence 
spectrum. The possibility of characterizing the spec- 
trum by a function with some arbitrary factors will 
then certainly be useful. 


The spectrology of turbulence will serve in the study 
of the flow in wind tunnels and will be interesting in 
meteorology. It seems to the author that its principal 
use may be the study of turbulence in the boundary 
layer or in wakes and especially in research concerning 
the mechanism of transition. 
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